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ABSTRACT 


Spins  of  sixteen  excited  states  of  30P  have  been 
investigated  in  conjunction  with  multipole  mixing  and 
branching  ratios  of  a  number  of  transitions  via  the 
reaction  28Si(3He,  p)30P  (Q  =  6.34  MeV) .  Spin  assign¬ 
ments  are  consistent  with  those  of  other  investigators 
in  all  cases  except  for  the  4.92  MeV  level.  The  branching 
ratios  for  this  level  are  also  in  disagreement  with 
those  previously  reported.  Assuming,  as  suggested  by 
the  present  results,  a  transition  to  the  0.68  MeV  level 
rather  than  to  the  0.71  MeV  level,  a  unique  assignment 
of  J  =  1  is  made  to  the  4.92  MeV  level.  A  probable  new 
assignment  of  J  =  1  has  been  made  to  the  3.74  MeV  level. 
Also,  the  previously  unstudied  2.84,  3.84,  3.93  and  4.73 
MeV  levels  have  been  assigned  possible  spin  values 
together  with  branching  ratios  and  multipole  mixing 
ratios  for  their  decay  transitions.  Branching  ratios 
for  transitions  from  the  3.74,  5.84,  3.93  and  4.50  MeV 
levels  are  in  serious  disagreement  with  those  of  previous 
studies.  Using  shell  model  wave  functions,  calculations 
of  lifetimes,  multipole  mixing  ratios  and  branching 
ratios  corresponding  to  five  low  lying  states  showed 
agreement  with  experimental  values  for  only  the  first 
two  excited  states. 
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CHAPTER  I  INTRODUCTION 

а .  Summary  of  Previous  work 

This  report  describes  a  study  of  level  spins  and 
gamma-ray  mixing  and  branching  ratios  for  levels  of 
30P  below  5  MeV  using  the  particle-gamma  angular 
correlation  method  of  Litherland  and  Ferguson  (Li  61) . 

The  first  major  investigations  of  the  self - conj ugate 
nucleus  30P  were  carried  out  using  the  32S(d,  a)30P 
reaction  (En  58)  and  the  29Si(p,  y)30P  reaction  (Va  58). 
Using  the  former  reaction,  Endt  and  Paris  accurately 
located  the  energies  of  30  excited  states  in  30P 
(to  ±10  keV)  in  a  region  of  excitation  up  to  5.8  MeV. 

They  also  observed  a  number  of  T  =  1  states  in  violation 
of  the  isospin  selection  rule  which  allowTs  only  T  =  0 
states  to  be  excited  in  this  reaction.  Using  the 
29Si(p,  y)30P  reaction,  Van  der  Leun  and  Endt  (En  58) 
made  the  following  spin,  parity  and  isospin  assignments 
'to  states  of  excitation  energy  E  (denoted  by  E  (J  ,  T)): 

Ground  (1+,  0);  0.68  (0+,  1);  0.71  (1+,  0);  1.45  (2+,  0) ; 
2.94  (2+ ,  1);  4.18  (2+,  1);  5.88  (2',  1);  5.96  (l',0); 

б. 23  Cl” ,  1) ;  and  6.27  (3" ,  0) . 

In  addition,  they  found 

the  isospin  selection  rule  forbidding  El  transitions 
between  levels  of  the  same  T  in  sel f- conjugate  nuclei 
to  be  valid  in  four  resonance  level  decays  that  went 
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predominantly  to  either  T  =  0  or  T  =  1  states.  Work 
on  the  same  reaction  by  Baart  et_  aJL.  (Ba  62)  gave  the 
following  spin  and  parity  assignments: 

Ground  (1  +  )  ;  0.71  (1  +  )  ;  1.45  (2  +  )  ;  1.97  (3+)  ;  2.54  (2+,  3+)  ; 
2.72  (2+ ,  (1+) )  ;  2.94  (2+) ;  3.02  (1+,  Cl'));  4.18  (2+,  C2")) 
5.89  (2');  6.00  (1*,  (1  +  ));  6.27  (2  +  )  ;  6.31  (3',  (3+)); 

6.52  Cl+,  (1’));  and  6.53  (2+) . 

The  less  probable  of  two  possible  assignments  has  been 
enclosed  in  parentheses.  Also  using  the  29Si(p,  y)30P 
reaction,  Moore  (Mo  63)  provided  additional  branching 
ratio  information  and  the  following  five  spin  and 
parity  assignments: 

2.54  (3  +  )  ;  4.23  (2+)  ;  6.84  (1  +  )  ;  7.04  (2  +  )  ;  and  7.21  (2  +  )  . 

% 

Further  work  by  Harris  and  Hyder  (Ha  66B,  Ha  67B)  on  the 
same  reaction  led  to  the  following  eight  unique  spin  and 
parity  assignments: 

0‘.  71  (1+)  ;  1.45  (2+)  ;  1.97  (3  +  )  ;  2.54  (3  +  )  ;  4.23  (4‘); 

4.47  (0+);  4.63  (3_);  4.93  (3‘);  and  7.04  (4‘). 

\ 

It  is  seen  that  the  spins  of  the  4.23  and  7.04  MeV  levels 
are  in  disagreement  with  the  assignments  by  Moore  (Mo  63) 
who  assumed  that  formation  of  such  a  strong  resonance 
by  g -  or  f - wave  protons  as  required  for  a  resonance  spin 
of  4  is  unlikely.  Using  a  2  c.c,  GeCLi)  detector, 
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Harris  and  Hyder  (Ha  67A)  investigated  decay  schemes  of 
7  resonances  in  the  29Si(p,  y)30P  reaction  and  provided 
an  accurate  energy  determination  (±2  keV  to  ±4  keV)  of 
17  bound  levels  of  30P.  They  found  a  weak  3%  transition 
from  the  J  =  4  ,  4.23  MeV  level  to  the  J77  =  1+  ground 
state  and  a  strong  transition  from  the^J77  =  3",  4.92 
MeV  level  to  the  J77  =  1+,  0.71  MeV  level.  Finally, 
lifetime  measurements  of  the  0.680  and  0.709  MeV  states 
by  Kennedy  et_  a^L.  (Ke  67)  yielded  values  of 
t  =  (1.55  ±  0.30)  x  10‘13  sec  and  x  £  3.0  x  10"12  sec, 
respectively.  They  also  reported  slightly  revised 
energies  of  677  ±  1  and  709.1  ±  1  keV  for  these  levels. 

b .  Model  Applications 

To  date  relatively  few  attempts  have  been  made  to 
interpret  30P  theoretically.  It  is  of  interest,  therefore, 
to  investigate  the  applicability  of  different  models 
by  considering  past  calculations  on  nuclei  in  this  mass 
region.  Successful  theoretical  interpretation  of  many 
nuclei  in  the  region  19  £  A  £  29  in  terms  of  a  rotational 
collective  model  has  been  possible.  For  example,  all 
observed  states  in  the  mirror  pair  Al25  and  Mg25  (Li  58) 
up  to  4  MeV  can  be  described  by  the  latter  model. 

Bromley  (Br  57)  showed  that  many  of 

the  properties  of  29Si  and  3.1P  can  be  explained  by  a 
rotational  collective  model.  Baart  et_  a_l.  (Ba  62) 
investigated  the  applicability  of  a  rotational  model 
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to  30P  and  found  it  inadequate  to  explain  the  energy 
levels  of  30P.  Picard  and  Pinho  (Pi  66)  applied  a 
strong-coupling  unified  model  with  a  Soper  two-body 
force  to  30P  but  satisfactorily  predicted  the  position 
of  only  the  first  two  excited  states.  Thankappan  and 
Pandya  found  that  a  weak- coupling  unified  model  explained 
electromagnetic  transitions  in  31P  better  than  did  a 
strong- coupling  unified  model  (Th  60).  Encouraged  by  this 
discovery  he  applied  the  former  model  to  30Si  and 
obtained  fair  agreement  with  the  observed  properties 
(Th  62)  .  Arima  (Ar  58),  on  the  basis  of  a  j-j  coupling 
shell  model  which  assumed  an  inert  28Si  core  and 
neglected  configuration  mixing,  investigated  nuclei  in 
the  mass  region  A  =  29  to  40.  In  the  case  of  3 °P  his 
predicted  energy  level  positions  showed  only  fair 
agreement  with  the  experimentally  observed  level  positions. 
Later  shell  model  calculations  by  Glaudemans  ejt  a_l.  (G1  64) 
which  included  configuration  mixing  produced  better 
agreement  with  observed  level  positions  in  the  case  of 
30P.  A  comparison  of  the  latter  model  with  observed 
radiative  properties  is  made  in  Chapter  IV.  Employing 
an  intermediate  coupling  model,  Bouten  e_t  a_l.  (Bo  67) 
calculated  energy  levels  as  a  function  of  a  parameter 
measuring  the  relative  strength  of  spin-orbit  and  central 
forces  for  almost  all  nuclei  in  the  mass  region  18  <  A  <  38. 
In  the  case  of  30P  their  predicted  energy  level  positions 
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show  only  fair  agreement  with  the  observed  energy  level 
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CHAPTER  II  EXPERIMENTAL  METHOD 

a .  Theory 

The  measurement  technique  used  in  this  work  involved 
the  detection  of  charged  particle  reaction  products 
emitted  near  180°  in  an  axially  mounted  annular,  silicon 
surface-barrier  detector.  Gamma  rays  in  coincidence 
with  these  protons  were  detected  by  aNal(Tl)  counter 
between  angles  of  0°  and  90°.  For  this  arrangement  the 
gamma- ray  angular  distributions  can  be  expressed  in 
terms  of  an  expansion  of  even  Legendre  polynomials  and 
are,  therefore,  symmetric  about  90°.  Because  of  the 
charged-particle  detector's  axial  symmetry,  the  angular 
distribution  will  also  be  symmetric  about  0°.  Litherland 
and  Ferguson  (Li  61)  have  shown  that  detection  of  the 
outgoing  particles,  y,  from  a  reaction  A(x.,  y)B*  by 
a  point  counter  on  the  beam  axis  limits  the  maximum 
number  of  magnetic  substates  of  the  residual  nucleus 
B*  that  can  be  populated  according  to 


lal  <  S ,  +  S  +  S 
1  'max  —  A  x  y 


where  a  is  a  magnetic  quantum  number  of  B*  and  S^, 
and  S  are  the  maximum  z-components  of  spin  of  the 

y 

target,  bombarding  and  outgoing  particles,  respectively. 
This  restriction  arises  basically  from  the  fact  that  the 
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orbital  angular  momenta  contained  in  plane  waves 
have  only  zero  projection  along  their  direction  of 
propagation.  A  further  condition  on  the  population 
parameters  results  from  the  alignment  of  the  state  B*. 

If  P (a )  denotes  the  population  of  the  a  substate 
relative  to  the  other  substates  then  alignment  as  opposed 

to  polarization  of  B*  implies  that  P(a)  =  P(-a).  Thus, 

for  the  reaction  28Si(3He,  p)30P,  (Q  -  6.344  MeV)  studied  in 

this  work  the  highest  substate  populated  is  given  by  |a|max  =  1. 

Consequently,  the  populations  of  the  three  substates 

a  =  -1,  0  and  1  are  fully  specified  by  the  ratio 

P(±1)/P(0)  together  with  the  normalization  condition 

£  P(a)  =1.  An  important  feature  of  the  above  technique 
a 

is  that,  generally,  B*  will  be  a  bound  state  of  definite 
spin  and  parity.  All  information  about  the  nature  of 
the  reaction  mechanism  is  contained  in  the • population 
parameters  which  are  treated  as  unknowns  in  the  analysis. 

The  effect  of  the  non-zero  size  of  the  particle 
detectors  is  to  populate  substates  higher' thanl a  I 

max . 

Litherland  and  Ferguson  investigated  this  problem  (Li  61) 
and  found  that  the  corrections  can  be  calculated  only 
with  a  knowledge  of  the  details  of  the  reaction.  They 
showed,  moreover,  that  corrections  in  lowest  order  are 
proportional  to  £2,  £  being  the  half-angle  subtended 
by  the  counter.  For  an  annular  detector  at  180°,  E,2  is 


. 
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replaced  by  (Z^  ~  where  Z^  and  ^  are  the  angles 

(measured  from  180°)  which  are  subtended  at  the  target 
by  the  inner  and  outer  edges  of  the  counter.  Since  the 
corrections  occur  because  of  the  nonzero  projection  of 
the  outgoing  orbital  angular  momentum  on  the  beam  axis 
it  is  expected  that  higher  angular  momenta  for  the  out¬ 
going  particles  will  lead  to  larger  corrections  to  the 
population  parameters. 

Following  the  formalism  of  Poletti  and  Warburton 
(Po  65)  the  angular  distribution  of  gamma-rays  arising 
from  the  decay  of  a  state  faf  to  a  state  ’bT  is  given  by 

W(e)  =  l  a  P  (cos  6)  =  l  pk(a)Fk(ab)QkPk(cos  9) 
k  k 

k  even  (1-1) 

where  6  is  the  angle  between  the  direction  of  emission 

of  the  gamma  rays  and  the  axis  of  alignment.  The  Qk  are 

attenuation  coefficients  for  the  gamma-ray  detector, 

f 

the  Pp(a)  are  statistical  tensors  which  describe  the 
alignment  of  the  initial  state,  and  the  Fk(ab)  depend 
specifically  on  the  gamma-ray  cascade  and  are  independent 
of  nuclear  alignment. 

The  statistical  tensors  pk(a)  can  be  expressed  in 
terms  of  a  weighted  sum  of  the  population  parameters 
of  the  2a  +  1  magnetic  substates  of  ’a’  by 

PkU)  =  I  Pk(a,  °opc°o 

a 


(1-2) 


'■ 
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Alignment  of  the  state  ’a’  implies  P(a)  =  P(-a),  and 

the  normalization  is  such  that  £  P(a)  =1.  If  the  sum 

a 

in  equation  2  is  limited  to  positive  a  then 


(a,  a)  =  (2 


,) 


(aaa-a 

ko) 

(aaa-a 

oo) 

(1-3) 


where  the  symbol  (aaa-a |ko)  represents  a  Clebsh-Gordan 
coefficient.  Values  of  Pp(a>  a)  for  integer  values  of 
a  between  1  and  6  and  for  half-integer  values  of  a 
between  3/2  and  11/2  are  listed  in  the  paper  of  Poletti 


and  Warburton  (Po  65)  .  The  F^.(ab)  are  given  by 


•  L  ’  -L  +  7T  -  7T 
l 


XLXL 


Fk(LL'ba)/j;  xL 


(1-4) 


where  L,  L*  are  integers  ranging  from  a-b  to  a+b  (with 
the  restriction  L,  L*  /  0),  tt  and  ir ?  are  associated  with 
L  and  L',  respectively;  tt  is  0  for  electric  and  1  for 
magnetic  radiation.  The  sign  convention  is  that  of 
Litherland  and  Ferguson  (Li  61).  The  quantities  and 
xT  f  are  given  by  ratios  of  reduced  matrix  elements  as 

L 


<b 

L  | 

a> 

<b 

L' 

a> 

<b 

L 

m 

|  a>  ’ 

XL  ’  <b 

L 

m 

a> 

The  reduced  matrix  elements  are  real  and  is  the  lowest 
allowed  value  of  L.  The  F^LL’ba)  are  given  by 


. 
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Fk(LL'ba)  =  O)13"3'1  [  (2L  +  1)  (2L'  +  1)  (2a+l)  ]^(L1L' -1  |ko) 

x  W(aa  LL' ;  ko)  (1-6) 

They  have  the  symmetry  property  F-^(LL’ba)  =  F^(L'Lba), 
and  the  normalization  Fo(LL’ba)  =  6  t .  The  triangular 
conditions  on  the  Racah  coefficient,  W(aa  LL ' ;  ko) , 
in  equation  6  limit  k  to  k  <  min(2L,  2L',  2a).  Values  of 
the  F^.(LL’ba)  are  tabulated  in  the  paper  of  Poletti 
and  Warburton  (Po  65) ,  and  in  the  Argonne  National 
Laboratory  report,  ANL  5324,  by  M.  Ferentz  and  N. 
Rozenzweig . 

Often  the  second  gamma  ray  in  a  cascade  a  b  -*  c 
provides  useful  information  either  in  conjunction  with  the 
first  transition  or  by  itself.  The  angular  distribution 
of  the  second  gamma  ray  is  given  by 


where 


W(e)  =  l  pk(a)Uk(ab)Fk(bc)QkPk(cos  0)  (1-7) 

•  k 


Uk(ab)  =  l  x£Uk (Lab) / £  x£  (1-8) 

L  L 


The  Uk(Lab)  are  given  by 


tt  ? T  W(abab;  Lk) 

Uk(Lab)  =  wjabab ;  Lo)' 


(1-9) 


so  that  Uo(Lab)  =  1  and  Uk(Lab)  =  Uk(Lba).  Equation  7 
can  be  generalized  to  give  the  angular  distribution  of 
the  n^^1  ‘gamma  ray  in  a  cascade  as 
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w(6)  =  I  Pk(a)Uk(ab)Uk(bc) . . .Uk(xy)Fk(yz)QkPk(cos  0) 

♦  k 

(1-10) 

which  contains  the  product  of  n-1  Uk  coefficients. 

Rarely  is  it  necessary  to  include  more  than  two 
multipoles  in  equations  4  and  8.  For  only  two  multipoles 
present  equations  4  and  8  become 

F,  (LLba)  -  (-)°2xF,(LL'ba)+x2F,  (L'L'ba) 

Fk(ab)  =  - - - £ -  (1-11) 

1  +  x2 

and 

U,  (Lab)  +  x2lh  (L ’ ab) 

Uk(ab)  =  - - -  (1-12) 

k  1  +  x2 

where  L  is  the  lowest  allowed  multipole  and  L'  =  L  +  1. 

The  denominator  1  +  x  provides  a  convenient  normalization. 

It  is  sometimes  useful  to  express  equations  11  and  12  in 

terms  of  an  angle  cf>  given  by  4>  =  arc  tan  x.  In  these 

equations  the  mixing  parameter  x  is  given  by 


<b 

L  +  l 

a> 

<b 

L| 

a> 

(1-13) 


so  that  the  x^  and  x^ ,  of  equation  5  have  become  x^  =  1 , 
xL,  =  x.  The  Litherland-Ferguson  sign  convention  takes 
a=0  for  an  EL+l/ML  mixture  and  a  =  1  for  an  ML+l/EL 
mixture.  The  Rose  and  Brink  convention  (Ro  67)  corresponds 
with  the  case  a  =  0  independent  of  the  type  of  multipole 
mixing.  The  latter  convention  is  used  throughout  this 


- 
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report.  Values  of  U^.(Lab)  are  tabulated  in  the  paper 
of  Poletti  and  Warburton  (Po  65) . 


b .  Method  of  Analysis 


The  general  problem  of  analysis  involves  fitting  the 
equation 

w(0)  =  I  Pk(a>  a)P(a)Fk(ab)QkPk(cos  6)  (1-1) 
ka 

to  a  measured  angular  distribution.  The  computer  program 
used  in  the  present  work  treats  the  population  parameters 
P(a)  as  unknowns  and.,  for  each  of  a  series  of  discrete 
values  of  the  mixing  ratio  (contained  in  Fk(ab))  and 
assumed  values  for  the  spin(s),does  a  least  -  squares 
fit  to  the  measured  points  in  the  angular  distribution. 

The  best  fit  corresponds  to  those  values  of  the  population 
parameters  which  give  the  lowest  value  of  x2  where 


i  (Y(e.)  -  WCe.))2 

n  l - V  (FT"— 

1  v 


(1-14) 


The  quantity  E(e^)  is  the  estimated  uncertainty  in  the 
measured  gamma-ray  yield  Y(0^)  at  angle  0^  and  N  is  the 
number  of  degrees  of  freedom.  One  condition  that  the 
function  defined  by  1-14  must  satisfy  for  it  to  be  a 
X2 -distribution  function  is  that  the  denominator  should 
equal  the  variance.  A  good  approximation  to  the  variance 
is  made  if  one  assumes  Poisson  statistics  and  takes  the 
number  of  counts  in  a  peak  for  E2(0^).  However ,  in  the 
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present  work  E2(6^)  includes  not  only  the  uncertainty 
in  the  peak  area  but  also  uncertainties  due  to  back¬ 
ground.  A  problem  arises  with  equation  14,  namely, 
that  of  defining  the  correct  number  of  degrees  of 
freedom.  For  example,  one  may  carry  out  the  fitting 
of  W C 8 )  to  the  observed  gamma-ray  yields,  Y(e^),  in  two 
different  ways.  One  way  would  be  to  vary  both  the 
population  parameters  and  the  mixing  ratio  together. 
Another  way  would  be  to  step  through  all  possible  values 
of  the  mixing  ratio  and  for  each  value  fit  W(e)  to  the 
yields,  Y(e^).  The  net  result  of  both  procedures  is 
a  solution  for  the  mixing  ratio  and  the  population 
parameters  on  the  basis  of  a  minimum  value  of  x2. 
However,  if  one  takes  the  number  of  degrees  of  freedom 
as  the  number  of  data  points  minus  the  number  of  free 
parameters  then  in  the  first  case  there  is  one  less 
degree  of  freedom  than  in  the  second  case.  Moreover, 
no  account  has  been  taken  of  the  effectiveness  in 
changing  W(e)  of  each  of  these  free  parameters.  For 
instance,  in  fitting  the  angular  distribution  from  a 
spin  zero  state  only  the  mean  of  the  data  points  is 
unknown.  In  fitting  a  nearly  isotropic  angular  distribu¬ 
tion  with  an  assumed  spin  of  1  for  the  initial  state, 
one  will  obtain  a  small  k  =  2  term  that  can  be  fitted 
with  a  number  of  different  mixing  ratios  and  population 


pa 


rameters.  This  effect  occurs  because  the  k  =  2  term  is 
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of  the  form 

A2  =  2F  (x° ,  x1,  ...  xp)[P(l)  -  P(0)] 

where  x  is  the  mixing  ratio  and  p  is  an  integer  such 
that  p  <_  3 .  It  can  be  seen  that  for  a  given  normalization 
(i.e.  a  given  value  of  P(0)  +  2P(1))  it  is  possible  to 
produce  a  small  A 2  for  a  wide  range  of  mixing  ratios. 

This  situation  is  equivalent  to  fitting  the  distribution 
with  the  Legendre  expansion  AQ  +  A2P2  (cos  0).  In  the 
latter  case  both  Aq  and  A2  are  uniquely  predicted  and 
the  number  of  degrees  of  freedom  can  be  unambiguously 
shown  to  be  N-2  (N  is  the  number  of  data  points).  Thus, 
the  three  parameters  x,  P(0)and  P(l)  and  the  number  of 
data  points  result  in  the  same  number  of  degrees  of 
freedom  as  the  two  parameters  Aq  and  A2  and  the  number 
of  data  points.  In  the  opposite  extreme,  however, 
there  are  many  cases  where  a  fit  leads  to  unique  values 
of  x,  P(0)  and  P(l)  in  which  case  all  three  of  these 
parameters  are  ’fully  effective’.  Determination  of  the 
correct  number  of  degrees  of  freedom,  r,  can  be  the 
difference  between  a  given  solution  being  at  the  10% 
confidence  limit  for  r  =  1  or  at  the  50%  confidence 
limit  for  r  =  2. 

In  the  present  work  the  population  parameters  are 
treated  as  unknowns  and  a  least  -  squares  fit  is  carried 
out  for  a  discrete  number  of  assumed  mixing  ratios. 
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The  number  of  degrees  of  freedom  is  taken  as  the  number 
of  data  points  minus  the  number  of  allowed  population 
parameters.  Clearly,  of  the  possible  alternatives, 
the  above  choice  will  give  the  highest  values  of  x2  and, 
therefore,  in  the  sense  of  rejecting  possible  solutions, 
this  choice  will  be  the  most  conservative  one.  More¬ 
over,  spin  zero  solutions  and  those  arising  from 
isotropic  angular  distributions  for  spin  1  will  have 
fewer  degrees  of  freedom  than  they  should  have. 

For  a  correct  solution  the  value  of  x2  should  be 
near  the  expected  value  of  unity.  The  probability  that 
for  a  correct  solution  the  value  of  x2  obtained  will  be 
greater  than  a  given  value  can  be  read  from  probability 
tables.  In  the  present  work,  solutions  which  did  not 
dip  below  the  0.1%  confidence  limit  were  rejected. 

In  the  reaction  28Si(3He,  p)30P  (Q  =  6.34  MeV)  both 
P(0)  and  P(l)  were  allowed.  Effects  of  the  finite  size  of 
the  annular  charged-particle  detector  were  estimated  by 
including  a  fraction  of  the  next  highest  substate.  One 
fit  was  done  with  P(±2)  =  0.025  P(0)  and  another  with 
P(±2)  =  0.05  P(±l).  The  curve  having  the  greatest  effect 
on  the  solution  was  plotted  together  with  the  curve 
for  P(2)  =  0.  Limits  on  the  mixing  ratio  solution  were 
taken  at  the  extremes  of  these  curves  as  shown  by  the 

a  D.  R.  Whitney,  Elements  of  Mathematical  Statistics, 

(Henry  Holt  and  Company,  New  York,  1959). 


. 
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Arc  fan  (mixing  ratio) 


Figure  2-1  Hypothetical  x2-versus  arc  tan  (mixing 

ratio)  least-squares  fit. 
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arrows  in  Figure  2-1.  The  value  of  x2  at  which  limits 
on  the  mixing  ratio  solution  were  taken  were  one  standard 
deviation  above  the  501  limit  if  the  minimum  was  at  or 
below  this  level.  If  the  minimum  was  above  the  50%  limit, 
it  was  first  normalized  to  the  50%  limit  and  the  un¬ 
certainty  limits  taken  as  mentioned  above.  In  cases 
for  which  excessive  normalization  was  required,  limits 
were  taken  at  the  0.1%  confidence  limit  without  first 
normalizing.  A  more  detailed  discussion  of  the  fitting 
theory  is  given  in  Appendix  I. 

c .  Description  of  the  Experimental  Apparatus 

Approximately  0.2  microamperes  of  6.2  MeV  singly- 
charged  3  He  beam  from  the  University  of  Alberta  Van  de 

■L 

Graaff°  was  used  to  bombard  a  self-supporting  target  of 

natural  silicon.  Charged  particles  scattered  into  the 

angular  range  of  170°  to  173°  were  detected  by  an  annular 

c 

silicon  surface  barrier  detector  in  the  180°  geometry. 

Gamma  rays  were  observed  between  angles  of  30°  and  90°  by 
a  (7.6  cm  *  7.6  cm)  Nal(Tl)  detector^.  The  charged- 

b  Model  CN,  High  Voltage  Engineering,  Burlington, 
Massachusetts,  U.S.A. 

c  Nuclear  Diodes,  10,000Q-cm,  750  microns  depletion  depth 
at  200  volts  bias. 

d  Harshaw  "Integral  Line’  Assembly  type  12S12/3. 
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particle  beam  originating  from  the  radio  frequency  ion 
source  first  passed  through  a  terminal  analyzing  magnet 
which  sorted  out  the  undesired  mass  components  from  the 
beam.  The  particles  were  then  accelerated  and  bent 
through  90°  by  a  double  focusing  sector  magnet  of  radius 
28.5  in.  Four  feet  from  the  exit  port  of  the  analyzing 
magnet  the  beam  passed  through  a  pair  of  energy  stabilizing 
slits.  The  current  balance  on  these  slits  was  maintained 
by  a  feedback  circuit  controlling  the  impedance  of  the 
corona  loading  of  the  accelerator  terminal.  This  system 
limited  the  energy  spread  of  the  beam  at  the  target  to 
less  than  5  keV.  The  beam  then  passed  through  a 
carbon  trap,  switching  magnet  and  quadrupole  lens.  Since 
the  straight  through  port  was  used  only  a  small  switching 
magnetic  field  was  necessary  to  steer  the  beam  into 
the  quadrupole  lens.  A  tantalum  collimator  with  a  1/4  in. 
diameter  hole  was  placed  at  the  entrance  to  the  quadrupole 
to  ensure  axial  alignment  with  the  beam.  Approximately 
9  feet  from  the  exit  of  the  quadrupole,  the  beam  passed 
through  a  second  tantalum  collimator  with  a  3/8  in.  diameter 
hole  backed  by  a  lh  in.  long  lead  cylinder  of  h  in. 
inner  diameter.  This  collimator  served  to  define  the  beam 
entrance  into  a  9  in.  long  tantalum  shielded,  copper 
carbon  trap  of  3/4  in.  inner  diameter.  This  carbon  trap 
minimized  carbon  buildup  on  the  collimators  and  target. 
Following  this  trap  were  four  more  tantalum  collimators 


* 
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holes  of  1/4,  1/8,  1/16  and  3/32  in.  in  that  order. 

The  first  three  collimators  were  each  backed  by  a  lh  in. 
long  lead  cylinder  of  1/4  in.  inner  diameter.  The  last 
collimator  was  placed  on  the  end  of  the  last  lead 
cylinder  and  functioned  as  an  antiscattering  shield. 

The  latter  four  collimators  were  spaced  over  a  length  of 
approximately  1  foot.  A  final  collimator  consisting  of 
a  tantalum  sheet  with  a  3/32  in.  diameter  hole  followed 
by  a  1  in.  long  lead  block  having  a  hole  slightly 
larger  than  that  in  the  tantalum  was  situated  in  the 
target  chamber  immediately  behind  the  annular  charged- 
particle  detector.  This  collimator  prevented  the 
detector  resolution  from  being  affected  by  intense  low 
energy  gamma  rays  that  would  otherwise  arise  from  the 
beam  hitting  the  tantalum  shield  on  the  back  of  the 
detector  and  spraying  on  the  sides  of  the  hole  through 
the  detector.  After  passing  through  the  target,  the 
beam  was  stopped  in  a  tantalum  lined  beam  catching  tube 
approximately  1  metre  long.  This  arrangement  reduced 
180°  elastic  scattering  into  the  charged-particle  detector 
to  a  negligible  amount. 

Angular  changes  in  the  gamma  detector  were  made  by 
moving  a  rotatable  arm  onto  which  the  Nal(Tl)  detector  was 
mounted.  The  arm  rotated  about  the  center  of  an 
aluminum  frame  which,  in  turn,  was  bolted  to  the 
concrete  floor.  The  entire  counter  support  assembly, 
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or  goniometer,  was  modeled  after  a  similar  design  at 
Brookhaven  National  Laboratory.  The  target  chamber  was 
supported  by  the  center  post  of  the  goniometer  and  was 
estimated  to  be  within  h  mm  of  the  axis  of  the  rotatable 
arm  that  held  the  gamma-ray  detector.  A  photograph 
of  the  target  chamber  mounted  on  the  goniometer  is 
shown  in  Figure  2-2. 

The  gamma-ray  detector  was  placed  at  a  distance 
of  10cm  from  the  target  and  had  a  resolution  of  6.2% 

(FWHM)  for  the  1.45  MeV  gamma  ray  of  30P.  Specifications 
for  this  detector  listed  a  maximum  gain  shift  of  less 
than  1%  for  a  counting  rate  change  from  1000  c/sec 
to  10,000c/sec  and  a  drift  of  less  than  1%  over  a  24 
hour  period  for  a  counting  rate  of  1000  c/sec.  An  1/8  in. 
lead  sheet  was  placed  over  the  detector  face  in  order 
to  reduce  the  intensity  low  energy  gamma  radiation. 

For  this  experiment  the  gamma  detector  was  set  at  angles 
of  31°,  45°,  60°  and  90°.  Interference  of  the  beam 
catching  tube  prevented  the  gamma  detector  from  being 
set  at  angles  less  than  31°.  Eight  runs  were  made 
with  two  runs  at  each  angle.  The  runs  were  typically 
6  hours  long. 

The  charged-particle  detector  was  placed  approximately 
2  cm  from  the  target.  Particles  were  detected  in  a 
sensitive  region  of  the  detector  having  a  1  mm  wide 
annulus  of  5.5  mm  inner  diameter.  For  10  MeV  protons 
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Figure  2-2  Photograph  of  the  Goniometer  and 

Target  Chamber 

Only  the  top  portion  of  the  goniometer  can  be  seen. 
The  final  lead  collimator,  annular  detector,  target  holder, 
part  of  the  beam  catching  tube  and  the  half-moon  shaped 
rotatable  arm  that  supports  the  gamma-ray  counter  can  all 
be  seen. 
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this  detector  exhibited  a  resolution  of  close  to  65  keV. 
This  figure  includes  the  effect  of  target  thickness. 

System  isotropy  was  checked  using  the  J77  =  0+, 

2.32  MeV  level  of  1  4N  formed  by  the  12C(3He,  p)14N 
(Q  -  4.779  MeV)  reaction.  The  difference  between  the 
normalized  number  of  counts  at  any  angle  and  the  mean, 
normalized  number  of  counts  was  zero  to  within  the 
2%  statistical  uncertainty. 

d  Electronics 

A  block  diagram  of  the  electronics  used  in  this 
experiment  is  shown  in  Figure  2-3.  Fast  coincidences 
between  the  Nal(Tl)  gamma- ray  and  annular  surface 
barrier  charged-particle  detectors  were  used  to  gate 
the  analog  pulses  from  these  detectors  into  two  independent 
analog- to-digital  converters  (ADC’s)  which  were  inter¬ 
faced  on-line  with  an  SDS  920  computer. 

Signals  from  the  gamma-ray  detector  were  pre- 

6  -f 

amplified  and  then  amplified  by  a  shaping  amplifier 
operating  in  the  double  delay  line  clipping  mode.  The 

bipol  ar  output  was  fanned  out  to  the  F  input  of  a 

4096  channel  ND-161F  Dual  Analog  to  Digital  Converter^  (ADC) 

through  an  adjustable  delay  line  for  pulse  height 
analysis,  and  to  a  cross-over  pickoff  circuit  for 
timing  purposes.  The  (20  nsec  wide,  0.5  volt  after  dif- 

e  OrteC  Model  113  Scintillation  Preamplifier, 

f  Ortec  Model  410  Multimode  Amplifier, 

g  Nuclear  Data  Inc.,  Palatine,  Illinois, 
h  Ortec  Model  407  Cross-over  Pickoff. 
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Figure  2-3 
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f erentiating)  output  of  the  latter  circuit  provided  the 

stop  input  to  a  time-to-amplitude  converter1  (TAC) . 

Signals  from  the  solid  state  charged-particle 

detector  were  preamplif ied^  and  then  fanned  out  to 
two  main  shaping  amplifiers  in  parallel.  One  of  these 
amplifiers,  operating  with  double  delay  line  shaping, 
provided  the  input  pulses  to  a  timing  single  channel 
analyzer  which  in  turn  fed  the  start  input  of  the 
TAC.  The  other  amplifier  was  operated  with  RC  shaping 
of  0.5  microseconds  integration  and  two  1  microsecond 
differentiations  to  present  a  suitable  pulse  shape  to 
the  M  input  of  the  ND-161F  ADC^  for  optimum  charged- 
particle  resolution. 

The  fast- coincidence  gating  logic  pulses  were 

obtained  from  the  TAC  by  setting  two  40  nanosecond  wide, 

contiguous  windows  around  its  output  by  means  of  a  dual 

0 

channel  analyzer  .  The  logic  output  from  this  unit 
corresponding  to  the  lower  window  corresponded  to  true 
plus  random  coincidences  and  was  used  to  gate  the  F  ADC 
and  interrupt  the  computer  for  sorting  of  the  coincident 
gamma- ray  and  charged-particle  pulses  into  certain 
regions  of  computer  memory.  The  logic  output  from  the 
upper  window  corresponded  to  random  coincidences  only 
and  was  used  to  interrupt  the  computer  for  sorting  of 
the  random  gamma- ray  and  charged-particle  pulses  into 

other  regions  of  computer  memory. _ __ 

i  Ortec  Model  437  Time- to-Pulse  Height  Converter, 

j  Ortec  Model  109  Transitor  Preamplifier 

lc  Ortec  Model  420  Timing  Single  Channel  Analyzer. 

F.  S.  Goulding  and  R.  A.  McNaught,  Nucl.  Instr. 

and  Methods  9(1960)  282. 
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The  gated  gamma  rays  associated  with  the  decay 
of  a  certain  level  or  excited  state  of  the  final  nucleus 
were  sorted  into  bins  by  the  computer  by  setting  a 
window  around  the  particle  group  feeding  that  level. 

Each  such  window  generated  a  bin  consisting  of  a  256 
channel  true-plus  -  random  gamma -ray  spectrum  followed 
by  a  256  channel  associated  random  gamma-ray  spectrum. 

Sixteen  such  bins  were  set  allowing  the  simultaneous 
investigation  of  the  decay  of  16  excited  states  of  the 
final  nucleus.  Also  was  recorded  a  256  channel  .  spectrum 
of  the  total  unsorted  true-plus-random  coincident  gamma 
rays,  followed  by  a  256  channel  spectrum  of  the  cor¬ 
responding  total  random  coincident  gamma  rays  and  the 
true-plus  -  random  coincident  charged-particle  spectrum 
(1024  channels)  followed  by  its  corresponding  random 
spectrum  (1024  channels) . 

Up  to  eight  bins  at  a  time  could  be  viewed  on  the 
computer  disp'lay  unit  while  the  experiment  was  in  progress. 

A  photograph  of  such  a  display  is  shown  in  Figure  2-4. 

The  data  for  each  run  were  stored  on  magnetic  tape 
and  processed  off-line  on  the  University  of  Alberta 
IBM  7040  Computer. 

Monitoring  of  the  experiment  was  accomplished  by 
recording  a  self-gated  charged  particle  spectrum  on  a 
separate  1024  channel  pulse  height  analyzer111  .  Normalization 
of  the  data  for  each  run  was  determined  from  the  area 

m  Technical  Measurements  Corp.  Model  210  256  Channel 

Pulse  Height  Analyzer. 
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Figure  2 

The 
channels 
the  last 


4  Photograph  of  8  Gamma-ray  Bins  Taken 
from  the  SDS  920  Computer  Display  Unit, 
trues  plus  randoms  are  in  the  first  256 
of  each  bin  followed  by  randoms  only  in 
256  channels. 
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of  selected  peaks  in  the  monitor  spectrum, 
e .  Target  Preparation 

The  self-supporting  targets  used  in  this  work 
were  prepared  by  evaporation  of  natural  silicon.  The 
procedure  involved  first  evaporating  NaCl  onto  glass 
slides.  These  slides  were  then  placed  at  varying 
distances  from  the  crucible  of  an  electron  gun  evaporator0. 
Silicon  was  evaporated  using  an  electron  current  of  100  ma 
at  3500  volts  in  a  vacuum  of  less  than  10~^  torr 
pressure.  A  slow  evaporation  was  necessary  to  avoid 
sputtering  from  the  molten  silicon  metal.  After 
evaporation ,  the  silicon  films  were  floated  off  the 
glass  slides  onto  distilled  water  and  then  picked  up 
on  aluminum  target  frames.  In  this  way  self-supporting 
targets  ranging  in  thickness  from  approximately  30  to 
200  yg/cm2  were  made.  Measurements  of  target  contamination 
and  thickness  were  obtained  by  elastically  scattering 
3  MeV  alpha  particles  from  several  of  the  targets.  For 
the  target  used  resolutions  of  56  and  47  keV  (FWHM) 
were  obtained  for  scattering  from  silicon  and  oxygen, 
respectively.  Using  an  estimated  value  of  35  keV  for 

o  Model  no.  980-0001  Single  Crucible  e-Gun  Evaporation 
source,  Varian  Associates,  Palo  Alto,  California. 
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detector  plus  electronics  resolution,  these  figures 
correspond  to  thicknesses  of  approximately  55  yg/cm2  and 
30  yg/cm2  for  silicon  and  oxygen,  respectively. 

All  other  contaminants  below  approximately  A  =  40  gave 
elastic  scattering  peak  areas  less  than  0.1%  of  the  area 
of  the  oxygen  peak. 

f .  Spectrum  Analysis 

Gamma  rays  can  interact  with  a  detector  crystal 
by  Compton  scattering,  photoelectric  absorption  or  pair 
production.  Ihe  probability  of  interaction  by  each  of 
these  processes  depends  on  the  gamma-ray  energy.  At 
low  energies  photoelectric  absorption  is  the  dominant 
process  and  gives  rise  to  a  single  peak  or  photopeak 
corresponding  to  absorption  of  the  full  gamma-ray  energy. 
Compton  scattering  is  present  at  all  energies  and  results 
in  a  broad  distribution  beginning  at  slightly  below  the 
full  energy  p'ealc  and  extending  down  to  zero  energy. 

It  is  also  possible  for  the  Compton  scattered  photon 
to  be  absorbed  in  the  crystal  and,  in  this  way,  contribute 
to  the  full  energy  peak.  The  relative  importance  of  the 
latter  process  depends  on  gamma-ray  energy  and  crystal 
size.  At  gamma-ray  energies  above  1.02  MeV  pair  produc¬ 
tion  becomes  energetically  possible.  The  probability 
for  interaction  by  pair  production  increases  with 
increasing  energy  while  that  for  photoelectric  absorption 

The  probability  for  interaction  by  Compton 


decreases . 
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scattering  also  decreases  with  increasing  energy  but 
less  rapidly  than  it  does  for  photoelectric  absorption. 

At  energies  above  approximately  1.5  MeV  three  other 
peaks,  in  addition  to  the  photopeak,  are  observed. 

These  peaks  result  from  the  interaction  of  the  photons 
by  pair  production  in  either  the  crystal  or  the 
surrounding  material.  The  positrons  resulting  from  the 
pair-production  process  annihilate  with  electrons 
producing  two  0.511  MeV  photons.  If  the  positrons 
are  produced  in  the  crystal,  one  or  both  of  the  annihila¬ 
tion  photons  can  escape  from  the  crystal.  If  only  one 
photon  escapes,  a  single  peak  0.511  MeV  lower  in  energy  than 
the  photopeak  is  produced  while  if  both  photons  escape, 
a  peak  1.02  MeV  lower  in  energy  than  the  photopeak  is 

t 

produced.  The  0.511  MeV  peak  results  from  the  detection 
of  annihilation  radiation  produced  in  the  surrounding 
material  and  not  from  the  primary  interactions  in  the 
detector.  Thus,  the  response  to  or  ’line  shape’  of  a 
given  gamma  ray  is  due  to  the  crystal  dimensions  and  a 
combination  of  the  above  three  basic  processes. 

To  obtain  branching  ratios  from  a  gamma-ray  spectrum 
one  must  extract  the  relative  intensities  of  the  different 
energy  gamma  rays  present.  With  spectra  obtained  by  Nal(Tl) 
detectors  these  intensities  are  usually  obtained  from 
photopeak  areas  corrected  for  photopeak  efficiency  as 


« 


it 
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a  function  of  energy^3.  Thus,  one  is  usually  concerned 
with  unfolding  contributions  from  a  mixture  of  a  number 
of  different  gamma -ray  'line  shapes'.  In  this  work 
gamma  rays  of  interest  extended  from  0.68  MeV  up  to  4.50 
MeV.  Standard  line  shapes  were  constructed  from  the 
gamma  rays  originating  from  the  1.45,  1.97,  2.72,  3.73, 

4.14  and  4..  93  MeV  levels  of  30P.  From  these  'line  shapes' 
curves  were  constructed  giving  the  ratio  of  the 
photopeak  height  to  the  heights  of  the  first  valley, 

Compton  shoulder,  first  escape  peak,  second  valley, 

second  escape  peak  and  asymptotic  Compton  tail,  respectively, 

as  a  function  of  gamma-ray  energy.  A  graph  of  full 

width  at  ha]f  maximum  (FWHM)  as  a  function  of  energy 

was  also  constructed.  Using  these,  self-consistent 

stripping  of  the  spectra  was  accomplished. 

It  was  found  that  slight  gain  changes  often 
occurred  between  runs,  partly  because  of  the  large 
count  rate  change  that  resulted  from  turning  off  the  beam 
and  partly  because  of  the  1/2  to  1  hour  of  shutdown 
necessary  to  carry  out  routine  operations.'  For  this 
reason  a  program  was  written  in  SDS  920  machine  language 
to  shift  the  counts  in  channel  x  of  a  spectrum  to  a  new 
position  x'  given  by  x'=Gx+B.  The  shifting  was  done 
so  as  to  make  corresponding  peaks  in  the  different  sets  of 
spectra  coincide.  This  program  uses  linear  interpolation  and 

p  R.  L.  Heath  et  al.,  IDO-17017  (1965)  A.E.C.  Research 


and  Development  Report. 
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was  found  to  preserve  areas  to  better  than  1  part  in  500.  A 
more  detailed  discussion  of  the  program's  format  and  general 
features  is  given  in  an  internal  report. g 

In  determining  branching  ratios  from  spectra  of 
gamma  rays  that  de-excite  one  level  only,  usually  only 
a  few  gamma  rays  are  present..  In  these  cases  the  analysis 
is  usually  relatively  simple.  Moreover,  a  self-consistent 
check  is  usually  available  in  that  corrected  intensities 
of  gamma  rays  in  a  given  cascade  must  be  consistent 
with  each  other  and  with  the  branching  ratios  of  lower 
lying  levels.  Thus,  analysis  must  start  from  the  lower 
excited  states  and  proceed  successively  towards  the  higher 
lying  levels. 

Photopeak  areas  in  this  work  were  determined  from 
plotted  spectra  by  summing  the  number  of  counts  under 
a  given  peak  and  subtracting  off  the  estimated  back¬ 
ground  due  to  higher  energy  gamma  rays.  The  line 
shapes  of  these  higher  energy  peaks  were  constructed 
using  the  parameters  of  standard  line  shapes  mentioned 
previously.  The  consistency  of  this  method  was  checked 
by  comparison  with  areas  of  several  peaks  as  determined 
by  a  modified  gaussian  peak  fitting  program  .  Uncertainties 


g  C.  Vermette,  University  of  Alberta  Internal  Report,  1967 

h  M.  Putnam  et  al.,  IDO-17016  (1965)  A.E.C.  Research 


and  Development  Report. 


. 


35- 


were  estimated  by  talcing  the  square  root  of  the  total 
number  of  counts  under  a  peak  including  background. 

An  additional  uncertainty  associated  with  defining 
the  background  was  also  included.  Uncertainties  in  the 
branching  ratios  were  computed  from  corrected  photopeak 
areas  according  to  the  expression 


d  B . 

sB±  =  II  — 1  c«Ap 

\l  i  J 


(1-B.)2(<5A.)2  +  B?(  y  6  A . ) 
l  v  l  i4 

1  J 
jVi 


(1-16) 


where  B^  is  the  branching  ratio  of  the  i*^  transition 
(J  B.  =  1),  A.  is  the  photopeak  area  corrected  for 

j  J  1 

relative  efficiency  of  the  gamma  ray  corresponding  to 
this  transition  and  SA.  is  the  uncertainty  in  A.. 
Branching  ratios  were  determined  from  the  summed  spectra 
of  runs  at  45°  and  60°.  Because  the  intensity  of  each 
gamma  ray  has  an  angular  distribution  of  the  form 


W(e)  =  Aq  +  A2P2(cos  6)  +  A4P4(cos  6)  (1-17) 

where  the  P^  (k  even)  are  Legendre  polynomials,  branching 
ratios  determined  in  the  above  mentioned  way  must  be 
suitably  corrected.  The  total  intensity  of  gamma-ray 
emission  into  a  solid  angle  of  4ir  steradians,  because 
of  the  axial  symmetry  and  symmetry  about  90°,  is  proportional 
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i  a  /  w(e)de  =  aq 

0° 

By  normalizing  both  the  data  and  W(e)  to  Aq  =  100 
correction  factors  to  the  branching  ratios  were  determined 
by  summing  the  values  of  the  data  at  45°  and  60°, 
dividing  by  2  and  taking  the  reciprocal.  Except  for  a 
few  distributions  having  large  A^  terms,  correction  factors 
were  usually  within  5%  of  unity. 


. 
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CHAPTER  III  EXPERIMENTAL  RESULTS 

The  following  is  a  detailed  discussion  of  experimental 
results.  Results  discussed  in  subsections  a  to  r  are 
summarized  in  tables  3-1,  3-2  and  3-3.  A  summary  of 
existing  information  on  levels  of  30P  below  5  MeV 
together  with  spin  assignments  and  branching  ratios 
obtained  in  the  present  work  is  given  in  figure  3-1. 

a .  Yield  Curves 

To  determine  the  best  3He  bombarding  energy  a 
charged-particle  yield  curve  using  doubly- charged  3He 
was  run  over  an  energy  range  of  4.7  to  10  MeV  in  steps 
of  100  keV.  The  greatest  variation  in  proton  yield 
took  place  in  the  neighborhood  of  5.5  to  6.5  MeV. 

Yield  curves  for  protons  corresponding  to  30P  levels 
of  the  most  interest  are  given  in  figure  3-2a  for  an 
energy  range  of  4.7  to  7.2  MeV.  Because  only  general 
trends  in  the  yield  curve  rather  than  accurately  known 
areas  were  required,  the  points  in  these  curves  taken 
as  the  maximum  number  of  counts  in  a  peak  measured 
from  the  background  level  for  10  microcoulombs 
of  integrated  beam  current.  A  bombarding  energy 
of  6.2  MeV  was  chosen  since  most  of  the  curves  in 
figure  3-2a  exhibit  a  maximum  or  near  maximum  at  this 
energy.  A  coincident  charged-particle  spectrum  taken  at 
a  3He  energy  of  6.2  MeV  is  shown  in  figure  3T2b. 
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Figure  3-1  Summary  of  Spin  Assignments  and  Branching 

Ratios  for  30P  Levels. 
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Table  3-1  Branching  Ratios  of  30P  States. 


Initial 

State 

Final 

State 

E 

Y  . 

(MeV) 

Branching  Ratios 
Present  Harris  Baart 

et  al.  et  al. 

Others 

1.453 

0 

1.453 

9  5±4 

>96 

100 

0.680 

0.773 

<  3 

0.709 

0.744 

5±4 

<  3 

1.974 

0 

1.974 

41±5 

34+2 

40 

56  (Va 

58) 

0.709 

1.263 

59  +  5 

53+2 

47 

44  (Va 

58) 

1.453 

0.521 

<4 

13±2 

13 

2.536 

0 

2.536 

95±4 

18  5 

95 

100  (Mo 

63) 

0.680 

1.8561 

,  Q 

c 

0.709 

1.827  J 

<  y 

o 

1.453 

1.083 

5  ±4 

1.974 

0.562 

<10 

2.722 

0 

2.722 

100 

18  0 

0.680 

2.042l 

<  7 

0.709 

2.013J 

^  Li 

1.453 

1.269 

<3 

1.974 

0.748 

<2 

2.839 

0 

2.839 

2  2  ±3 

100 

0.709 

2.130 

50+4 

1.453 

1.386 

28±  4 

2.939 

0 

2.939 

10+6 

19±2 

15 

0.680 

•2.259 

43+6 

37±2 

42 

1.453 

1.486 

4  7+6 

44±2 

43 

3.020 

0 

3.020 

<1 

<  5 

11 

0.680 

2.340 

100 

>_95 

89 

3.736 

0 

3.736 

46±6 

18  0 

53  (Mo 

63) 

0.680 

3.056 

29±5 

1.453 

2.283 

18±6 

4  7  (Mo 

63) 

2.939 

0.797 

7±3 

3.836 

0 

3.836 

6±3 

100  (Mo 

63) 

0.680 

3.156 

18±9 

1.453 

2.383 

16±4 

2.939 

0.897 

60±6 

• 

' 
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Initial 

State 

3.928 

4.142 

4.182 

4.231 

4.501 

4.624 

4.734 

4.921 


Table  3-1  (cont’d) 


Final  E  Branching  Ratios 

Present  Harris  Baart  Others 

State  (MeV)  et  al.  et  al . 


0 

0 

<6 

0.680 

0.709 

3.24  8' 
3.219 

10±6 

1.453 

2.475 

35±  6 

(100) 

2.939 

0.989 

55±6 

0 

4.142 

70±4 

88±4 

87 

(Mo 

63) 

0.709 

3.433 

14±5 

12  +  4 

1.453 

2.689 

16±5 

12 

(Mo 

63) 

0 

4.182 

19±3 

6 

0.709 

3.473 

81±3 

62 

2.939 

1.243 

32 

0 

4.231 

<10 

3  ±1 

1.453 

2.778 

12±  5 

1.974 

2.257 

69±4 

69±  2 

64 

(Mo 

63) 

2.536 

1.695 

19±4 

2  8±  2 

26 

(Mo 

63) 

0 

4.501 

42±5 

70  +  8 

66 

0.709 

3.792 

<9 

1.453 

3.048 

58±5 

30+8 

3.020 

1.481 

34 

0.680) 

*3.915 

<11 

14 

(Mo 

63) 

0. 709J 

(Mo 

63) 

1.453 

3.171 

5  9  ±  7 

71  +  3 

86 

1.974 

2.650 

14  +  6 

5±  2 

2.939 

1.685 

2  7±5 

24±3 

0 

4.734 

9±5 

0.680 

4.054' 

|  1 5±  5 

0.709 

4.025 

1.453 

3.281 

18±5 

2.939 

1.794 

58  +  6 

0.680 

4.241 

56±7 

0.709 

4.212 

63±  5 

1.453 

3.468 

37±5 

2.939 

1.982 

19±3 

4.231 

0.690 

2  5  ±4 

■ 


. 

' 

«• 
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Table  3-2  A 2  and  Coefficients  from  Legendre 

Polynomial  Fits  to  the  Angular  Distribution. 


Energy 

Level 

(MeV) 

Decay 

Trans ition 

A2 

: 

> 

X2 

1.453 

1.453+0 

0.83±0.03 

0.0010.05 

0.7 

1.974 

1.974+0 

0.29+0.10 

-0 .3810. 16 

1.7 

1.974+0.709 

0.29±0.08 

-0.1510.14 

0.2 

0.709+0 

0.37+0.07 

0.10+0.12 

0.7 

2.536 

2.536+0 

0.29+0.10 

-0.52+0.15 

0.2 

2.722 

2.722+0 

0.42+0.10 

0.36+0.14 

<0.1 

2.839 

2.839+0 

0.4710.11 

-0.10+0.17 

<0.1 

2.839+0.709 

0.50+0.15 

-0.37+0.23 

0.4 

3.020 

3.020+0.680 

0.74+0.14 

-0.1110.23 

0 . 5 

0.680+0 

0.01+0.03 

-0.0010.05 

0.2 

3.736 

3.736+0 

1.16+0.22 

-0.22+0.31 

0.5 

0.680+0 

-0. 0310. 06 

-0.01+0.11 

0.5 

3.836 

3.836+2.939 

0.3210.12 

-0.2310.19 

0.1 

0 . 680+0  ■ 

0.10±0.10 

0.05+0.15 

1.4 

3.928 

3.928+2.939 

0.39+0.09 

-0.0610.16 

0.4 

0.680  n 

0.  709-"0 

-0.0710.08 

-0.10+0.15 

<0.1 

4.142 

4.142+0 

-0.43+0.03 

-0.07+0.05 

1.0 

4.501 

4.501+0 

-0.0H0.03 

0 . 09  +  0 . 05 

2.1 

4.501+1.453 

0.0410.13 

-0. 1010'.  21 

<0.1 

1.453+0 

-0.15+0.07 

-0.06+0.12 

1.8 

4.624 

4.624+1.453 

0.27+0.07 

-0.05+0.11 

0.5 

4.624+2.939 

0.10+0.27 

0.42+0.43 

<0.1 

4.734 

4.734+2.939 

-0.27+0.07 

0.1410.11 

0.6 

2.939+0.680 

+0.2510.17 

-0.26+0.25 

0.1 

4.921 

4.921+2.939 

-0.11+0.11 

-0.09+0.17 

0.1 

4.921+0.680 

-0.55+0.03 

0.0110.06 

<0.1 

0  •  680  n 

0.690  0 

-0.03+0.04 

-0.01+0.06 

0.3 

I 
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Table  3-3  Summary  of  Spin  Assignments  and  Mixing 

Ratio  Solutions. 


Energy 

Decay 

Assumed 

X2 

Level 

Transition 

Spin 

min 

Solutions  of  Mixing  Ratio, x 

(MeV) 

Sequence 

1.453 

1. 453+0 

1+1 

.7 

-2.36<x<-.45  or  .58<x<1.73 

2+1 

.7 

. 14<x< 1.19 

1.974 

1.974+0 

2  +  1 

1.7 

x> 6 . 3  or  x<- 0 . 40 

3+1 

1.7 

xl“  *  +0  17 

x<-6.3  or  x>19  or  x=- (0 . 49_ q * ^ 

1.974+0.709 

2+1 

.3 

3+1 

.2 

x=+C0.2i:°0-32°6)  or  x=+(2.5*3‘2) 

1 . 974+0 . 709' 
0 . 7  0  9->  0 

■  3+1+1 

.8 

x(0.709^0)  =  -(0.25^-2^)  or 

=  _  (4  oi  +  19  ) 

1  *  -  2 . 2  8 J 

1.974+0 

j 

2.536 

2.536+0 

2+1 

.2 

x<-  .  42  or  x> 28 

3+1 

.2 

r  ot  +  0.39s  _a  +  1.97^ 

x=+  (  .  23  _  q  2Q)  or  x=+(2.36_.7q  ) 

2.722 

2.722+0 

1+1 

3.5 

(x=Undef ined) 

,  +°° 

2+1 

.1 

x="(n  -6.3^ 

+  3  3  N 

1.5 

x=" 34- .10^ 

2.839 

2.839+0 

1+1 

.4 

+  n  q? 

xf_0  or  x=+  (1 . 00 _  q  ’  ^g)  or  xl23 

2+1 

.1 

x=-(5.7_3  ^Jor  x=-(0.65_q 

• 

3+1 

.1 

+ 00  +  D  ^  1 

x=+  (5 . 7  _  2  6)or  x=-(0.03_0*15)or 

xf-  5 .  7 

2.839+0.709 

1+1 

2.4 

x£.  04 

3.2 

x=+ (i . 00_  68  ) 

2+1 

.4 

-231x1-0.67 

3+1 

.4 

x=+C0.03^‘Jg)or  x=+(4.5^56) 

3.020 

3.020+0.680 

1+0 

.6 

x  =  0 

0.680+0 

0+1 

.2 

x  =  0 

3.736 

3.736+0 

1+1 

.5 

-8  lxl-0.13 

2+1 

2.5 

-4.331x1-0.60 

3+1 

2.2 

x<  .  04  or  xl4 . 5 

3.836 

3.836+2.939 

1+2 

.9 

xl~.21  or  .231x11.54 
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Energy 

Level 

(MeV) 


3.928 


4.142 


4.501 


4.624 


Table  3-3  (cont’d) 


Decay 

Assumed 

x2 

Transition 

Spin 

Sequence 

min 

Solutions  of  Mixing  Ratio, 

2+2 

.1 

x~C2.9:^) 

.9 

-  0 . 58<x<  +  0 . 38 

3+2 

1.1 

x=-(0.47*°;^)  or 

2.0 

x--  (3 • 7_2>3) 

4+2 

.1 

x<-19  or  x=+ (0 . 11±0 . 23) 

r  ,  _  +  OO 

or  x=+(6.3_^  g) 

3.928+2.939 

1+2 

.4 

x<-  .  12  or  x>5.7 

2  +  2 

-  .4 

0 . 55<x<6 . 3 

3+2 

.4 

x=(0._>06) 

1.2 

x=+(3. 49^o [ 74^ 

4.142+0 

1+1 

2.2 

-29  <x<-.05  or  .31<x<3.17 

2+1 

1.0 

x=+t2-i4-o:3o) 

1.6 

x=+  C • 1 5  _ . 22) 

4.501+0 

0+1 

2.9 

x=0 

1+1 

3.5 

(x=Undef ined) 

2+1 

2.3 

,  f  r  7  +3.84>. 
x  +^5*7  - 1 .34^ 

3.8 

x=-  (.181.11) 

4.501+1.453 

0  +  2 

.1 

XE  0 

1+2 

.2 

(x=Undef ined) 

2  +  2 

.1 

x=+(5-7  - 1 . 16} 

4.501+1.4531  „  -  , 
1.453+0 

.2 

1.4 

x=-(*19-!o8) 

(x (1 . 453+0) =Unde fined) 

1  +  2+1 

.7 

(x (1 . 453+0) =Unde fined) 

4.624+1.453 

1  +  2 

.6 

x< -  .05 

,  ^+.44. 

1.0 

x  =  +  (. 7 5 _  ^ 3Q) 

2+2 

.  6 

x=+  (  .  62 _  13)or  x  =  +  (5 . 7 _ 3  Q 

3+2 

.8 

x=- ( . 04± . 08) 

2.0 

x-+ (5 • 7_3  . Q) 

' 


f  ' 
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Energy 

Level 

(MeV) 


4.734 


4.734 


4.921 


Table  3-3  (cont’d) 

Assumed  "  x2 

Transition  Spin  min  Solutions  of  Mixing  Ratio, x 

Sequence 


4 . 624+2 . 939 

1+2 

LO 

• 

O 

(x=Undef ined) 

2+2 

.1 

x=-  (5- 7!4.9)°r 

x> -  0 . 2  9 

3+2 

.1 

x=- (0.16_  07  )°r 

.3  or  x<4 . 0 

4.734+2.939 

1+2 

1.3 

x<  0  or  x> 3 . 5 

3.0 

. 28<x< 1 . 57 

2+2 

.6 

. 65<x<  9 . 5 

3+2  ' 

1.3 

x=+I.02±0.05) 

1.0 

x=+(4-4-i;4) 

4.734+ 

2.939+0 

1+2+0 

0.7 

lxlf_0.78  or  |  x  | 

>1.9 1 

- 

X  ” 

2+2  +  0 

3+2+0 

0.1 

|x|<0.70 

;  x(4. 734+ 
2.939) 

0.5 

Mil- 54 

J 

4. 734+2.939] 

1+2  +  0 

.  7 

x=- C.27_>09) 

2. 939+0. 680  J 

2+2+0 

1.2 

x=4.33_3.68) 

3+2+0 

.  7 

For  x=0.02 

4 . 92 1+ . 709 

1+1 

.1 

- 4 . 5<x< .23  or  . 

34<x<  2 . 2  5 

0 

2  +  1 

.1 

x=+(. llt’11)  or 

x-+(l . 73+0.32) 

4.921+.680 

1+0 

.1 

x  =  0 

4.921+2.939 

1+2 

.2 

(x=Undef ined) 

2+2 

.3 

.  rl+.OS. 
x  ^  *  51- . 26 

.1 

|  x  |  >_1 4  . 
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Figure  3-2a  Yield  Curves  for 
The  points  correspond  to 
of  counts  in  a  peak  with  backg 
for  10  microcoulombs  of  charge 


10  States  of  30P. 
the  maximum  number 
round  subtracted 
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Figure  3-  2b  Coincident  charged-particle  spectrum 

from  °He  on  a  natural  silicon  target 
at  6.2  MeV  bombarding  energy. 
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b.  The  0.68  and  0.71  MeV  Levels 


Because  of  the  small  energy  separation  of  this 
doublet,  it  was  not  possible  to  study  these  levels  by 
observation  of  the  gamma  rays  in  coincidence  with  the 
protons  forming  these  states.  However,  correlations 
of  gamma  rays  de-exciting  these  levels  could  be  studied 
by  observing  known  0.68  and  0.71  MeV  gamma  rays  arising 
from  the  decay  of  higher  excited  states.  In  this  way, 
the  0.68  ->  0  MeV  transition  resulting  from  the  3.02  ■> 
0.68  ->  0  MeV  cascade  was  found  to  be  isotropic  with  a 
X2  value  of  0.2.  Furthermore,  in  a  Legendre  Polynomial 
fit  suspected  0.68  MeV  gamma-rays  arising  from  the  decay 
of  the  3.74,  3.84  and  4.92  MeV  levels  all  gave  values  of 
zero  for  and  to  within  a  standard  deviation. 

From  the  large  value  of  A^  obtained  in  the  Legendre 
Polynomial  fit  to  the  0.71  ->  0  MeV  transition  that 
arises  from  the  1.97  ->  0.71  0  MeV  cascade,  it  was 

possible  to  rule  out  spin  zero  for  the  0.71.  MeV  level. 

In  figure  3-3  curve  A  shows  a  simultaneous  fit  to  the 
1.97  ->  0.71,  0.71  ->  0  and  1.97  ->  0  MeV  transitions 
which  yields  mixing  ratio  solutions  of  x  =  -  CO. 25  +  q*  21^ 

+  1  Q 

or  x  =  -(4.01  £  28^  f°r  ^  MeV  transition 

while  curve  B  shows  a  fit  to  the  0.709  -*  0  MeV 
transition  alone  for  x(1.97  ±  0.71)  =  0.21.  In  this 
fit,  values  of  the  mixing  ratios  for  the  1.97  ->  0.71 
and  1.97  +  0  MeV  transitions  were  obtained  from 
separate  fits  to  these  transitions.  Although  x2_fits 
to  the  1.97  -0.71  MeV  transition  gave  a  large  region 
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Figure  3-3  Angular  Distribution  and  x2-Versus  Arc 

tan  x  for  the  0.71  -*  0  MeV  Transition 
that  Arises  from  the  1.97  ->  0.71  ->  0 


MeV  Cascade. 
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of  acceptable  values  for  the  mixing  ratio  of  this 
transition,  only  the  value  x  =  0.21  corresponding  to  the 

minimum  with  the  largest  quadrupole  component  was  used. 

Our  results  are  in  fair  agreement  with  the  values  x  = 

-(0.48  ±  0.04)  and  x  =  -(2.0  ±  2)  measured  by  Harris 
and  Hyder  (Ha  67A) . 

c .  The  1.45  MeV  Level 

The  gamma-ray  spectrum  corresponding  to  this  level 
(figure  3-4)  reveals  a  small  transition  to  the  0.7 
MeV  doublet.  From  the  accurately  known  position  of 
the  0.68  and  0.71  MeV  gammas  as  obtained  from  fits  to 
the  0.7  MeV  gammas  of  other  level  decays,  it  was 
possible  to  accurately  fix  the  position  where  a  0.77 
MeV  peak  would  occur.  The  weakness  of  any  such  gamma 
ray  indicated  that  the  1.45  MeV  level  decays  mainly 
to  the  0.71  MeV  level,  rather  than  to  the  0.68  MeV 
level.  The  branching  ratios  obtained  are  (95  ±  4)%  to 
ground  and  (5  ±  4)1  to  the  0.71  MeV  level.  These 
values  include  estimates  of  overlap  with  the  tail  of 
the  proton  peak  corresponding  to  the  0.7  MeV  doublet. 

Least- squares  fits  shown  in  figure  3-5  failed  to 
distinguish  between  J  =  1  and  J  =  2  assignments. 

Previous  assignments  (Va  58,  Ba  62,  Ha  66A)  to  this 
level  are  all  J  =  2+.  The  mixing  ratio  solution  for 
J  =  2,  0.14  <  x  £  1.19,  is  consistent  with  the  value 
measured  by  Harris  and  Hyder  (Ha  66b)  of  x  =  +(0.21  ±  0.03). 
Uncertainties  have  been  taken  at  the  17%  confidence  limit. 


d.  The  1.97  MeV  Level 
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Figure  3-4  Gamma- ray  Spectrum  for  the  1.45  MeV 

Level  of  30P. 

The  points  in  the  simultaneous  randoms  spectrum 
are  indicated  by  open  circles. 
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Figure  3-5  Angular  Distribution  and  x2_Versus  Arc 

tan  x  for  the  1.45  ->  0  MeV  Transition. 
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29Si(p,  y)30P  reaction,  both  Moore  (Mo  63)  and  Harris 
and  Hyder  (Ha  67A)  have  found  a  13%  branch  to  the  1.45 
MeV  level.  Baart  et_  al_.  assigned  an  uncertain  14% 
branch  to  the  1.45  MeV  level  from  analysis  of  the  695 
keV  resonance  decay  scheme  while  Van  der  Leun  and  Endt 
found  no  branch  to  this  level  after  examination  of  the 
same  resonance.  The  gamma-ray  spectrum  for  the  1.97 
MeV  level  shown  in  figure  3-6  reveals  a  small  1.45 
MeV  peak  that  is  attributed  to  a  3%  overlap  with  the 
tail  of  the  proton  peak  corresponding  to  the  1.45  MeV 
level.  The  branching  ratio  results  of  (41  ±  5)%  and 
(59  ±  5)%  for  transitions  to  ground  and  the  0.71  MeV 
level,  respectively,  are  in  poor  agreement  with  the 
values  (34  ±  2)%,  (53  ±  2)%  and  (13  ±  2)%  for  branching 
to  the  ground  0.71  and  1.45  MeV  states,  respectively, 
measured  by  Harris  and  Hyder  (Ha  67A) .  Least- squares 
fits  to  the  1.97  MeV  transition  shown  in  figure  3-7 

f 

excluded  a  J  =  1  assignment  at  the  1%  confidence  limit. 
Shown  in  figure  3-8  is  the  least-squares  fit  to  the 
1.97  ->  0.71  MeV  transition  which  gave  possible  solutions 
for  J  =  1,  2  and  3.  Previous  investigators  (Ba  62, 

Ha  66A)  have  assigned  J  =  3  to  this  level.  For  the 
1.97  +  0.71  MeV  transition  a  J  =  3  least-squares  fit 
gave  mixing  ratios  of  x  -  +(0.21_q  ?^)  and  x  -  +  (2.5_-^2 
with  an  83%  probability  while  a  J  =  3  fit  to  the  1.97  -> 
MeV  transition  gave  an  allowed  region  for  the  mixing 
ratio  of  x  >  -0.14  with  a  99.9%  probability.  No  other 
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Figure  3 

The 
are  indi 


-6  Gamma- ray  Spectrum  for  the  1.97  MeV 

Level  of  30P. 

points  in  the  simultaneous  randoms  spectrum 
cated  by  open  circles. 
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Figure  3-7 


Angular  Distribution  and  x?~Versus  Arc 
tan  x  for  the  1.97  ->  0  MeV  Transition. 
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Figure  3-8  Angular  Distribution  and  x2-Versus  Arc 

tan  for  the  1.97  •»  0.71  MeV  Transition. 
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mixing  ratios  for  these  transitions  are  available  for 
comparison  at  the  present  time. 

e .  The  2.54  MeV  Level 

As  seen  in  figure  3-2b,  the  proton  peak  corresponding 
to  the  2.54  MeV  level  reveals  the  presence  of  a 
contaminant.  A  3.1  MeV  gamma  ray  in  the  spectrum  of 
figure  3-9  indicates  that  the  reaction  of  3He  with  this 
contaminant  has  a  higher  Q  value  than  that  for  the 
28Si(3He,  p)30P  reaction.  This  fact  together  with  the 
large  unexplained  1.26  MeV  gamma  ray  in  the  spectrum 
of  figure  3-10  suggests  the  contaminant  is  2 9Si (2 9Si (3He ,  p)31P 
Q  =  10.18  MeV).  Calculations,  however,  show  that  the 
particle  peak  would  correspond  to  a  state  at  6.76 
MeV  in  31P.  Such  a  state  has  not  yet  been  reported. 

Other  unexplained  gamma  rays  occur  at  0.92  MeV  and 

probably  0.7  MeV.  Analysis  of  the  other  gamma-ray 

0 

peaks  yielded  branchings  of  (95  ±  4 ) %  and  (5  ±  4)% 
to  the  ground  and  first  excited  states,  respectively. 

These  results  are  in  good  agreement  with  previous 
measurements  (Ha  67A) .  The  weak  intensity  of  a  possible 
1.83  MeV  gamma  ray  indicates  a  branching  of  less  than 
5%  to  the  0.7  MeV  level.  A  leas t - squares  fit  shown  in 
figure  3-10  excluded  J  =  1  with  a  99.7%  probability.  Both 
J  =  2  and  J  =  3  gave  acceptable  fits  yielding  mixing 
ratio  solutions  of  x  <-0.42,  x  >  28  and  x  =  +(0.23_q  ^0^’ 
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Figure  3-9  Gamma- ray  Spectrum  for  the  2.54  MeV 

Level  of  30P. 

The  points  in  the  simultaneous  randoms  spectrum 
are  indicated  by  open  circles. 
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Figure  3-10  Angular  Distribution  and  x2-Versus  Arc 

tan  x  for  the  2.54  0  MeV  Transition. 
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+  1.97 

x  =  +(2.36_0*70)  for  J  =  2  and  3,  respectively.  For 
J  =  4,  pure  octupole  radiation  is  rejected  at  the  0.1% 
confidence  limit.  The  previous  assignment  for  this 
level  as  measured  by  Moore  (Mo  63)  and  Harris  and  Hyder 
(Ha  66B)  is  J  =  3+. 

f .  The  2.72  MeV  Level 

The  2.72  MeV  level  gamma- ray  spectrum  shown  in 
figure  3-11  indicates  a  branch  of  100%  to  ground. 

Previous  measurements  by  Harris  and  Hyder  (Ha  67A) 
indicate  a  80%  branching  to  ground  while  Baart  et_  al . 

(Ba  62)  obtained  a  100%  to  ground  branching.  Previous 
work  by  Baart  et  al_.  (Ba  62)  gave  this  level  a  J  ~  2+, 

(1  +  )  assignment.  x2_  fits  to  the  2.72  0  MeV 

transition  of  this  work  reject  a  J  =  1  assignment  with 
a  96%  probability  (see  figure  3-12).  Mixing  ratio 
solutions  for  J  =  2  are  x  =  -(11.4  an<i  x  ~  (0«  34_  jj '  ^q) 

with  the  former  value  being  4  times  more  probable  than 
the  latter  one.  In  view  of  the  positive  parity  in¬ 
dicated  by  the  previous  assignment  (Ba  62)  the  results 
of  this  work  suggest  a  predominantly  E2  transition. 

g  ‘  The  2.84  MeV  Level 

The  gamma  spectrum  corresponding  to  this  level  is 
shown  in  figure  3-13.  The  branching  ratios  are 
(22  ±  3)%,  (50  ±  4)%  and  (28  ±  4)%  to  the  ground, 
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Figure  3-11  Gamma-ray  Spectrum  for  the  2.72  MeV 

Level  of  30P. 

The  points  in  the  simultaneous 
are  indicated  by  open  circles. 
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Figure  3-12  Angular  Distribution  and  x2 -Versus  Arc 

tan  x  for  the  2.72  ->  0  MeV  Transition. 
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Figure  3-13  Gamma-ray  Spectrum  for  the  2.84  MeV 

Level  of  30P. 

The  points  in  the  simultaneous  randoms  spectrum 
are  indicated  by  open  circles. 
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0.7  MeV  and  1.45  MeV  levels,  respectively.  These 
values  are  in  disagreement  with  the  1001  to  ground 
branch  indicated  by  Moore  (Mo  63).  The  position 
of  the  0.7  MeV  gamma  peak  indicates  that  the  major 
part  of  the  transition  to  the  0.7  MeV  doublet  goes  to 
the  0.71  MeV  level.  Because  of  extensive  overlap 
with  the  2.94  MeV  level  in  the  proton  spectrum  extraction 
of  a  reliable  correlation  for  the  0.71  -*  0  MeV  transi¬ 
tion  was  not  possible.  However,  the  2.84  ->  0  MeV  and 
2.84  -*■  0.71  MeV  transitions  were  not  seriously  affected 
by  this  overlap  since  the  2.94  MeV  level  spectrum  does 
not  give  strong  gamma  rays  in  regions  near  the  2.13 
and  2.84  MeV  peaks.  Fits  to  the  2.84  -+  0.71  MeV 
transition  shown  in  figure  3-14  gave  minimum  values 
of  x2  below  the  0.11  limit  for  J  =  1,  2,  3.  However, 
the  solutions  for  J  =  2  and  3  are  3  times  more  probable 
than  those  for  J  =  1.  Fitting  both  the  2.84  0.71 

MeV  and  2.84  ->  0  MeV  transitions  simultaneously  failed 
to  provide  additional  discrimination.  In  the  J  =  2  fit 
to  the  2.84  ->  0.71  MeV  transition  a  mixing  ratio  of 
-23  <  x  £  -0.67  was  obtained  while  for  the  J  =  3 
fit,  two  possible  solutions  of  x  =  +  (0 . 03*q ' J 9)  and 

x  =  +(4.5  +  ^5.)  were  found.  The  2.84  ->  0  MeV  transition 

+  17 

shown  in  figure  3-15  yielded  solutions  of  x  =  -(5.7_3 >3) 

x  =  - C0.65to*i|D  for  J  =  2  and  x  =  +(5-7-2.6) 
x  =  - (0  .  03H  0  *  ?  1)  and  x<_-  5  .  7  for  J^3.  No  previous  spin  assignments 

“  U  •  1  J 
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4  Angular  Distribution  and  x 2 -Versus  Arc 
tan  x  for  the  2.84  +  0.71  MeV  Transition. 
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Figure  3-15  Angular  Distribution  and  x2-Versus  Arc 

tan  x  for  the  2.8  4  ->  0  MeV  Transition. 
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h .  2.94  MeV  Level 

Overlap  of  the  proton  peaks  corresponding  to  the 
3.02  and  2.84  MeV  levels,  and  the  nature  of  the  decay 
scheme  combined  to  make  it  impossible  to  extract 
reliable  correlations  for  the  gamma  rays  from  this 
level.  However,  branching  ratio  determinations  from 
the  gamma-ray  spectrum  of  figure  3-16  were  possible. 
Decays  to  the  ground,  the  0.68,  and  the  1.45  MeV 
levels  were  found  to  take  place  with  the  branching  of 
(10  ±  6)%,  (43  ±  6)%,  and  (47  ±  6)%,  respectively. 

These  values  are  in  fair  agreement  with  those  of 
previous  investigators  (Ba  62,  Mo  63,  Ha  67A) .  It  was 
not  possible  in  this  work  to  determine  to  which  of  the 
levels  of  the  doublet  at  0.7  MeV  the  2.94  MeV  level 
decayed.  Thus,  based  on  previous  work  (Va  58,  Ba  62, 

Mo  63,  Ha  67A)  a  transition  to  the  0.680  MeV  level  was 
assumed.  As  pointed  out  by  A.  R.  Poletti  it  is  sur¬ 
prising  that  an  E2  transition  to  the  0.660  MeV  level  is 
favoured  over  an  isotopic  spin  allowed  Ml  transition 
to  the  0.709  MeV  level.  On  the  basis  of  Weisskopf 
estimates  an  r(E2)/r(Ml)  ehnancement  of  ~104  is  required 
for  this  to  be  so.  Van  der  Leun  and  Endt  (Va  58) 
have  made  a  J71  =  2+,  T  =  1  assignment  to  the  2.94  MeV 
level.  Baart  et  aJL  (Ba  62)  have  confirmed  the  latter 
Jn  assignment. 
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i.  The  3.02  MeV  Level 

The  gamma-ray  spectrum  for  this  level  shown  in 
figure  3-17  indicates  a  20  to  30%  overlap  of  the  charged 
particles  feeding  this  level  with  those  feeding  the 
2.94  MeV  level.  This  value  is  supported  by  estimates 
based  on  the  coincident  charged-particle  spectrum. 
Fortunately  only  15%  of  the  total  area  of  the  strong 
2.34  MeV  gamma  peak  consisted  of  the  2.26  MeV  gamma 


s 


A.  R.  Poletti,  Private  communication,  1967. 
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Figure  3-16  Gamma- ray  Spectrum  for  the  2.94  MeV 

Level  of  30P. 

The  points  in  the  simultaneous  randoms  spectrum 
are  indicated  by  open  circles. 
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Figure  3-17  Gamma- ray  Spectrum  for  the  3.02  MeV  Level 

of  30P. 

The  points  in  the  simultaneous  randoms  spectrum 
are  indicated  by  open  circles. 
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ray  arising  from  the  decay  of  the  2.94  MeV  level  to  the 
0.68  MeV  level.  An  approximate  correlation  of  this 
2.26  MeV  gamma  ray  was  used  in  subtracting  off  this 
peak  from  the  2.34  MeV  peak  produced  by  the  3.02  +  0.68 
MeV  transition.  Least-squares  fits  to  the  angular 
distribution  of  the  3.02  0.68  MeV  transition  (see 

figure  3-18)  rejected  all  spins  except  J  =  1  at  the 
0.1%  confidence  limit.  The  value  of  0.6  for  x2  confirms 
a  previous  J  =  1  assignment  by  Baart  ejt  aJL .  (Ba  62)  . 

The  100%  branch  to  the  0.68  MeV  level  obtained  in 
the  present  work  is  in  good  agreement  with  the  95%  result 
of  Harris  and  Hyder  (Ha  67A)  but  shows  poorer  agreement 
with  the  89%  and  11%  branching  to  the  0.68  MeV  and 
ground  states,  respectively,  as  reported  by  Baart  et_  a_l. 

(Ba  62). 

j .  The  3.74  MeV  Level 

The  gamma  spectrum  of  figure  3-19  shows  a  (46  ±  6)% 
of  branch  to  ground,  (29  ±  5)%  to  the  0.68  MeV  level, 

(18  ±  6)%  to  the  1.45  MeV  level  and  (7  ±  3)%  to  the  2.94 
MeV  level.  Fits  to  the  crossover  transition  as  shown 
in  figure  3-20  give  J  =  1,  2  and  3  as  possible  solutions 
with  J  =  1  :  2  :  3  in  the  likelihood  ratio  of  5  :  1  .  5  :  1. 
The  mixing  ratio  for  the  J  =  1  solution  of  this  transition 
Is  in  the  region  -8  <  x  <  -0.13.  No  pievious  spin 
assignments  have  been  made  to  this  level. 
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Figure  3-18 


Angular  Distribution  for  the  3.02 
MeV  Level  of  30P. 
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Figure  3-19  Gamma-ray  Spectrum  for  the  3.74  MeV 

Level  of  30P. 

1  he  points  in  the  simultaneous  randoms  spectrum 
are  indicated  by  open  circles. 
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Figure  3-20  Angular  Distribution  and  x2“Versus  Arc 

tan  x  for  the  3.74  0  MeV  Transition. 
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k.  The  3.84  MeV  Level 

No  previous  information  exists  for  this  level. 

Analysis  of  the  gamma  spectrum  of  figure  3-21  yields 

a  branching  of  (60  ±  6)2  to  the  2.94  MeV  level, 

(16  ±  4)%  to  the  1.45  MeV  level,  (18  ±  9)%  to  the  0.7  MeV 
doublet  and  (6  ±  3)%  to  the  ground  state.  The  0.7  MeV 

peak  position  suggests  that  most  of  the  3.84  +  0.7 
MeV  transition  is  to  the  0.68  MeV  level.  Least- 
squares  fits  to  the  3.84  ->  2.94  MeV  transition  shown  in 
figure  3-22  reject  spin  zero  at  the  0.1%  confidence 
limit.  The  minimum  values  of  x2  for  J  =  1,  2,  3  and  4 
are  0.9,  0.1,  1.1  and  0.1  respectively.  An  upper  limit 
of  20  nsec  can  be  set  on  the  lifetime  of  this  level 
by  virtue  of  the  time  coincidence  conditions.  Longer 
lifetimes  would  give  appreciable  tail  on  the  upper 
side  of  the  time  coincident  peak  that  would  appear 
in  the  randoms  window  and,  thus,  lead  to  a  smaller  than 
normal  true  to  randoms  ratio.  An  estimate  of  the 
magnitude  of  this  effect  can  be  made  by  constructing 
a  time  coincidence  peak  for  a  state  with  a  20  nsec 
lifetime  from  an  assumed  peak  shape  for  a  state  with 
a  short  lifetime.  The  construction  shown  in  figure 
3-23  involves  drawing  curves  for  the  exponential 
decay  given  by 

N  (t )  =  Nq  exp  (-  Iq-)  (3-1) 


where  N(t)  is  the  number  of  nuclei  in  the  3.84  MeV 
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Figure  3-21  Gamma -ray  Spectrum  for  the  3.84  MeV 

Level  of  30P. 

The  points  in  the  simultaneous  randoms  spectrum 
are  indicated  by  open  circles. 
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state  at  time  t  (in  nanoseconds) .  At  equal  intervals 
on  the  time  axis  curves  defined  by  equation  3-1  were 
drawn  taking  the  values  of  No  to  be  the  number  of 
counts  at  each  of  the  intervals.  Next  the  number  of 
counts  at  a  given  interval  was  obtained  by  summing  the 
contributions  from  each  of  the  preceding  intervals. 

As  seen  in  figure  3-23  approximately  30%  of  the  amount 
of  the  coincidence  peak  in  both  windows  lies  in  the 
randoms  window.  The  shape  of  the  time  coincidence 
peak  for  a  state  with  a  short  lifetime , which  is  dependent 
mainly  on  electronics , may  be  slightly  different  than 
the  assumed  shape.  For  this  reason  a  conservative 
estimate  of  20%  rather  than  30%  is  used.  However, 
even  20%  of  the  time  coincidence  peak  in  the  randoms 
window  should  result  in  a  clearly  visible  nonrandom 
background  spectrum.  Because  of  the  20  nsec  lifetime 
limit  and  the  large  branch  to  the  0.68  MeV  spin  zero 

level,  a  J  =  4  assignment  requires  |M(E4)|2  >  1.5  x  103 
W  u  (r  =  1.2  fermis)  for  the  3.84  +  0.68  MeV  transition 
and  can,  therefore,  be  ruled  out.  In  the  latter  estimate 
two  standard  deviations  have  been  allowed  for  M3  in 
the  branching  ratio. 
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£.  The  5.93  MeV  Level 

This  level  has  been  observed  only  weakly  in  proton 
capture  work.  A  100%  branching  to  the  0.71  MeV  level 
as  suggested  by  Harris  and  Hyder  (Ha  67A)  is  clearly 
not  in  accord  with  the  observed  gamma  spectrum  as  seen 
in  figure  3-24.  Analysis  of  this  spectrum  indicates 
a  (55  ±  6)%  branch  to  the  2.94  MeV  level,  a  (35  ±  6)% 
transition  to  the  1.453  MeV  level  and  a  (10  ±  6)% 
transition  to  the  0.7  MeV  level.  Least  -  squares  fits 
to  the  3.93  ->  2.94  MeV  transition  shown  in  figure  3-25 
give  spin  possibilities  of  J  =  1,  2  or  3.  The  mixing 
ratio  solutions  corresponding  to  these  spin  values 
together  with  values  of  x2  for  the  best  fits  are 

9 

listed  in  table  3-1.  Because  of  the  difficulty  in 
extracting  reliable  correlations  for  the  other  transi¬ 
tions  no  other  peaks  in  this  spectrum  were  fitted. 

m.  The  4.14  and  4.18  MeV  Level 


From  the  charged-particle  spectrum  of  figure  3-2b 
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Figure  3-24  Gamma-ray  Spectrum  for  the  3.93  MeV 

Level  of  30P. 

The  points  in  the  simultaneous  randoms  spectrum 
are  indicated  by  open  circles. 
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Figure  3-25  Angular  Distribution  and  x2 -Versus  Arc 

tan  x  for  the  3.93  ->  2.94  MeV  Transition. 
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it  is  seen  that  the  4.14  -  4.18  doublet  and  the  4.23 
MeV- peaks  are  not  resolved.  it  was,  however,  possible 
to  distinguish  the  gamma  rays  of  the  4.14  -  4.18 
MeV  doublet  from  those  of  the  4.23  MeV  level  by  setting 
two  adjacent  bins  in  the  charged -part icle  spectrum  such 
that  the  4.14  -  4.18  MeV  doublet  occupied  most  of 
one  bin  and  the  4.23  MeV  peak  occupied  most  of  the 
other  bin.  The  resulting  gamma-ray  spectrum  for  the 
4.14  -  4.18  MeV  bin  is  shown  in  figure  3-26.  This 
procedure  enabled  branching  ratios  to  be  determined. 
From  the  charged-particle  spectrum  it  is  evident  that 
the  4.18  MeV  level  is  much  more  weakly  excited  than 
the  4.14  MeV  level.  By  subtracting  off  the  effect  of 
overlap  with  the  4.23  MeV  level,  branching  from  the 
4.14  MeV  level  of  (70  ±  4) % ,  (14  ±  5)%  and  (16  ±  5)1 
for  transitions  to  the  ground  state  the  0.71  MeV  and 
1.45  MeV  levels,  respectively,  were  obtained.  Harris 
and  Hyder  (Ha  67A)  report  a  branching  of  (88  ±  4)%  to 
ground  and  (12  ±  4) %  to  the  0.71  MeV  for  this  state 
For  the  decay  from  the  4.18  MeV  level  they  obtained 
(19  ±  3)%  and  (81  ±  3)%  for  transitions  to  ground  and 
the  0.71  MeV  levels,  respectively.  The  weak  (14  ±  5)% 
branch  to  the  0.71  MeV  level  indicated  by  our  results 
is  additional  evidence  that  the  4.18  MeV  level  is  very 
weakly  excited.  x2  fits  to  the  4.14  MeV  gamma-ray 
correlation  shown  in  figure  3-27  gave  J  =  1  and  2  as 
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Figure  3-26  Gamma- ray  Spectrum  for  the  4.14  -*•  4.18 

MeV  Levels  of  30P. 

The  points  in  the  simultaneous  randoms  spectrum 
are  indicated  by  open  circles. 
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Figure  3-27  Angular  Distribution  and  x2 -Versus  Arc 

tan  x  for  the  4.14  -*  0  MeV  Transition. 
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the  only  possible  solutions  with  J  =  2  being  5  times 
more  likely  than  J  =  1.  Two  possible  mixing  ratios 
are  obtained  for  J  =  2.  Taking  into  account  the  finite 
size  effect  of  the  annular  detector,  these  are  x  =  +  (0 . 15^  ‘  .j", ) 
to  the  0.1%  confidence  limit  and  x  =  +(2.14  q’^q)  to  the 
17%  confidence  limit.  The  latter  solution  is  about  1.9 
times  more  probable  than  the  former.  This  level  has 
previously  been  assigned  J  =  2  (Ha  66B)  . 

n.  The  4.23  MeV  Level 

This  level  has  previously  been  assigned  a  spin  of 
4”  with  a  (3  ±  1)%  branch  to  ground,  a  (69  ±  2)% 
transition  to  the  1.97  MeV  level,  and  a  (28  ±  2)% 
transition  to  the  2.54  MeV  level  (Ha  61k).  The  branching 
ratios  of  the  present  work  are  (19  ±  4)%  to  the  2.54 
MeV  level,  (69  ±  4)%  to  the  1.97  MeV  level  and  (12  ±  5)% 
to  the  1.45  MeV  level.  Figure  3-28  clearly  shows  a 
branching  to  the  1.45  MeV  level  not  found  in  other 
work  (Mo  63,  Ho  61k).  Unfortunately,  because  the  gamma 
rays  from  this  level  were  so  much  less  intense  than  those, 
of  the  4.14  MeV  level,  it  was  not  possible  to  obtain 
correlations  using  the  method  employed  for  the  4.14 
MeV  level. 

6.  The  4.50  MeV  Level 


Analysis  of  the  gamma-ray  spcctium  ox  xiguic  3  2./ 
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Figure  3-28  Gamma-ray  Spectrum  for  the  4.23  MeV 

Level  of  30P. 

The  points  in  the  simultaneous  randoms  spectrum 
are  indicated  by  open  circles. 
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Figure  3-29  Gamma- ray  Spectrum  for  the  4.50  MeV 

Level  of  30P. 

The  points  in  the  simultaneous  randoms  spectrum 
are  indicated  by  open  circles. 
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gave  branching  ratios  of  (42  ±  5)1  and  (58  ±  5)%  to 
the  ground  state  and  the  1.45  MeV  level,  respectively. 
The  low  intensity  of  the  0.7  MeV  gamma  ray  seen  in 
the  spectrum  of  figure  3-29  rules  out  the  possibility 
of  an  appreciable  4.50  3.02  MeV  transition  as  reported 

by  Baart  et_  al_.  (Ba  62)  .  Branching  ratios  from  the 
present  work  are  in  poor  agreement  with  those  of  Harris 
and  Hyder  (Ha  67A)  who  obtained  a  (70  ±  8 ) %  branch  to 
the  ground  state  and  (30  ±  8 ) %  branch  to  the  1.45  MeV 
level.  Fits  to  the  4.50  ->  0,  4.50  1.45  and  1.45  ->  0 

MeV  transitions  as  seen  in  figures  3-30,  3-31  and  3-32, 
respectively  gave  acceptable  solutions  for  J  =  0,  1, 
and  2.  To  obtain  better  discrimination,  a  simultaneous 
fit  of  the  1.45  ->  0  and  the  4.50  ->  1.45  MeV  transitions 
was  carried  out  by  stepping  through  possible  mixing 
ratios  for  the  4.50  ->  1.45  MeV  transition.  The 
minima  in  x2  for  these  fits  is  plotted  as  a  function 
of  mixing-ratio  solutions  for  the  4.50  1.45  MeV 

transition  in  figure  3-33.  In  the  resultant  fits 
minima  in  x2  were  obtained  at  a  mixing  ratio  for  the 
1.45  -*  0  MeV  transition  of  x  =  0.70  in  all  cases. 
However,  in  the  latter  fit  x2  values  for  all  values  of 
the  mixing  ratio  for  the  1.45  ->  0  MeV  transition  were 
below  the  17%  confidence  limit.  Thus,  the  mixing  ratio 
for  the  4.50  ->  1.45  MeV  transition  was  not  determined 
for  J  =  1.  The  x2  fit  for  J  =  1  and  x(4.50  ->  0)  =0.18 
is  shown  in  figure  3.34.  It  should  be  emphasized  that 
any  of  the  values  of  x(4.50  ->  0)  are  possible  so  that 


-117- 


Figure  3-30  Angular  Distribution  and  x2 -Versus  Arc 

tan  x  for  the  4.5  0  0  MeV  Transition. 
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Figure  3-31  Angular  Distribution  and  x2 -Versus  Arc 

tan  x  for  the  4.50  ->  1.45  MeV  Transition. 
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Figurc  3-32  Angular  Distribution  for  the  1.45  +  0 

MeV  Transition  Arising  from  the  4.50  + 
1.4  5  ->  0  MeV  Cascade. 
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Figure  3-33  Minimum  in  x 2 -Versus  Arc  tan  x  (Mixing 

Ratio  of  4.50  ->  1.45  MeV  Transition) 
for  1.45  0  MeV  Transition  that  Arises 

from  the  4.50  -*  1.45  -*  0  MeV  Cascade. 
J(4.50)  =  1  is  Assumed. 


the  minimum  in  x2  for  J  =  1  may  range  from  0.67  to  0.93. 

As  seen  from  figure  3-31,  for  J  =  2  the  two  possible 
mixing  ratio  solutions  are  x=+ (5 . 7+f  Jg)  or  x  =  -  (0 . 19*°  *  Jg)  . 
Simultaneous  fits  for  both  x(4.50  ->  1.45)  =  5.7  and 
x  (4 . 50  ->  1.45)  =  -0.19  were  made  to  the  4.50  +  1.45  and 
1.45  ->  0  MeV  transitions.  For  the  known  mixing  ratio 
of  the  1.45  -*  0  MeV  transition  of  x  =  0.21  ±  0.03  (Ha  66B)  ) 
the  fits  for  x(4.50  -t  1.45)  =  5.7  gave  a  minimum  value 
of  x2  of  -4.1  while  those  for  x(4.50  ->•  1.45)  =  —0.19  gave 
a  minimum  x2  value  of  3.0.  Only  the  former  fit  is 
shown  in  figure  3-34.  Thus,  a  spin  2  assignment  can  be 
rejected  at  the  1%  confidence  limit. 

It  has  been  proposed  by  Endt  and  Paris  (En  58) 
that  this  level  is  the  T  =  1  analogue  to  one  of  the 
doublet  levels  of  30S,  at  3.78  MeV.  In  addition, 

Harris  and  Hyder  (Ha  67A)  observed  that  the  4.50  MeV 
level  is  excited  by  a  primary  transition  from  the  J  =  2, 

1470  keV  resonance  in  the  29Si(p,  y)30P  reaction. 

Because  of  these  observations  together  with  the  fact 
that  the  4.50  MeV  level  decays  strongly  to  the  J  =  2, 

1.45  MeV  level,  Harris  and  Hyder  (Ha  67A)  have  proposed 
that  this  state  is  the  T  =  1,  J  =  1  analogue  to  the 
3.77  MeV  level  of  30Si.  Moreover,  they  have  also 
found  a  J  =  0+,  T  =  1,  4.468  MeV  level  that  is  probably 
the  analogue  to  the  3.79  MeV  level  of  30Si  (Ha  67B) .  Our 
results  support  a  J  =  1  assignment  to  the  4.50  MeV 


level . 
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Figure  3-34  Simultaneous  Least  -  Squares  Fit  to  the 

4.5  0  -*  1.4  5  and  1.45  0  MeV  Transitions 

for  the  4.50  MeV  Level  of  30P. 
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p.  The  4.62  MeV  Level 

Branching  ratios  of  (59  ±  7)%,  (14  ±  6)1  and  (27  ±  5 )% 
to  levels  at  1.45,  1.97  and  2.94  MeV  were  obtained  from 
the  spectrum  of  figure  3-35.  As  can  be  seen  in  figure 
3-36,  acceptable  x2  fits  to  the  4.62  ->  1.45  MeV  transitions 
were  obtained  for  J  =  1,  2  and  3.  A  J  =  0  assignment 
is  rejected  with  a  99.7%  probability.  A  J  =  3" 
assignment  has  been  made  to  this  level  by  Harris 
and  Hyder  (Ha  66B) .  For  J  =  3  mixing  ratios  of  x  = 

-(0.04  ±  0.08)  and  x  =  +(5.7*2*q)  in  a  probability 
ratio  of  4  :  1  were  obtained.  Based  on  a  J  =  3 
assignment  the  mixing  ratio  solutions  indicate  that  the 
4.62  -*  1.45  MeV  transition  is  probably  pure  dipole. 

Shown  in  figure  3-37  is  the  angular  distribution  of 
and  x2  fit  to  the  4.62  ->  2.94  MeV  transition.  Table 
3-3  summarizes  the  results. 

f 

q.  The  4.73  MeV  level 

The  branching  ratios  for  this  hitherto  unstudied 
level  obtained  from  the  gamma  spectrum  shown  in  figure 
3-38  are  (9  ±  5)%,  (15  ±  5)%,  (18  ±  5)%  and  (58  ±  6)% 
to  ground,  the  0.7  MeV  doublet,  1.45  and  2.94  MeV 
levels.  x2  fits  to  the  4.73  ->  2.94  MeV  transition 
shown  in  figure  3-39  gave  J  =  1,  2  and  3  as  possibilities. 
Simultaneous  least-squares  fits  to  the  4.73  -*  2.94  and 
2.94  +  0.68  MeV  transitions  gave  possible  assignments 
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Figure  3-35  Gamma-ray  Spectrum  for  the  4.62  MeV 

Level  of  30P. 

The  points  in  the  simultaneous  randoms 
spectrum  are  indicated  by  open  circles. 
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Figure  3-36  Angular  Distribution  and  x2 -Versus  Arc 

tan  x  for  the  4.6  2  ->  1.4  5  MeV  Level 

of  30P. 
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Figure  3-37  Angular  Distribution  and  X2-Versus  Arc 

tan  x  for  the  4.6  2  -*  2.94  MeV  Transition. 
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Figure  3-38  Gamma-ray  Spectrum  for  the  4.7 

Level  of  30P. 

The  points  in  the  simultaneous  randoms 
are  indicated  by  open  circles. 
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Figure  3-39  Angular  Distribution  and  x2 -Versus  Arc 

tan  x  for  the  4.73  -*  2.94  MeV  Level  of 

30P. 
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of  J  -  1,  2  and  3  in  the  probability  ratio  2:1:2 
(figure  3-40).  Mixing  ratio  solutions  for  these  possible 
assignments  are  given  in  table  3-3.  In  the  simultaneous 
fit  to  the  4.73  ->  2.94  -*  0.68  MeV  cascade  only  the 
solution  x  =  (0.02  ±  0.05)  obtained  from  the  fit  to 
the  4.73  2.94  MeV  transition  alone  was  used  since 

the  fit  to  the  2.94  0.68  MeV  transition  from  this 

cascade  excluded  the  other  solution  (see  fig.  3-40b). 

r .  The  4.92  MeV  Level 

The  gamma  spectrum  corresponding  to  this  level  is 
shown  in  figure  3-41.  Previous  work  by  Harris  and 
Hyder  (Ha  66B,  Ha  67A)  on  this  level  indicate  a  63% 
branching  to  the  0.71  MeV  level.  The  results  of  the 
present  work,  however,  suggest  a  strong  branching  to 
the  0.68  MeV  level  rather  than  to  the  .71  MeV  level. 

Peak  positions  of  the  .7  MeV  gamma  ray  arising  from 
the  decay  of  the  4.92  level  as  obtained  by  fits  with 
a  gaussian  at  the  eight  angles  at  which  data  were  taken 
agreed  within  0.10  channels  with  the  position  of  the 
known  0.68  MeV  gamma  ray  that  results  from  the  3.02  MeV 
level  decay.  x2  fits  t0  an  assumed  4.92  +  .71  MeV 
transition  shown  in  figure  3-43  yielded  possible  spin 
assignments  of  1  or  2.  This  result  is  in  disagreement 
with  the  J  =  3"  assignment  of  Harris  and  Hyder  (Ha  67A) . 
Least-squares  fits  to  the  4.92  ->  0.68  MeV  transitions 
yield  a  unique  J  =  1  assignment  for  the  4.92  MeV  level. 


Figure  3- 40a  Simultaneous  Leas t - Squares  Fit  to  the 

^  •  73  ->  2.94  and  2.94  ->  0.68  MeV  Transitions. 
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Figure  3-40b  Angular  distribution  and  x2 -Versus 

Arc  tan  x  for  the  2.94  0.68  MeV 

Transition  Arising  from  the  4.73  -* 
2.94  ->  0.68  Cascade. 
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Figure  3-41  Gamma-ray  Spectrum  for  the  4,92  MeV 

Level  of  30P. 

The  points  in  the  simultaneous  randoms  spectrum 
are  indicated  by  open  circles. 
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Figure  3-42 


t 


Angular  Distributions  of  the  4.92  ->  0.68 
and  4.92  ->  2.94  MeV  Transitions. 
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Figure  3-43  Least  -  Squares  Fit  to  an  Assumed  4.92  +  0.71 

MeV  Transition  for  the  4.92  MeV  Level  of 

30P. 
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Branching  ratios  of  (56  ±  7) % ,  (19  *  3)1  and  (25  ±  4)% 

to  the  0.68,  2.94  and  4.23  MeV  levels,  respectively, 

were  obtained.  A  surprising  feature  of  the  decay 

scheme  is  the  strong  branch  to  the  4.231  MeV  level 

which  has  been  assigned  J  =  4  by  Harris  and.  Hyder 
(Ha  66B)  .  Such  a  transition  requires  |M|2  >  7  *  103  for 

a  J  =  1  assignment  to  the  4.92  MeV  which  is  clearly 

incompatible  with  a  J  =  4  assignment  to  the  4.23  MeV 

level . 
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CHAPTER  IV  MODEL  CALCULATIONS 

Model  calculations  by  Glaudemans  et  al.  ( G 1  64) 
have  led  to  predicted  energies  and  configurations  for 
a  number  of  energy  levels  of  nuclei  in  the  mass  region 
between  29Si  and  40Ca.  Their  method  involved  treating 
15  matrix  elements  for  the  two-particle  potentials 
together  with  binding  energies  to  the  2 8 Si  core  of  a 

nucleon  in  the  2s^  and  ld^  shells  as  unknown  parameters. 

2  2 

Of  these  known  levels  the  states  of  30P  used  were  the 
0.68,  0.71,  2.54,  2.72  and  2.94  MeV  levels.  In  the 
present  work,  the  phase  consistent  formalism  of  Rose 
and  Brink  (Ro  67)  was  used  to  calculate  mixing  ratios 
and  transition  probabilities  of  some  of  the  states  of 

o  n 

P  predicted  by  Glaudemans  ert  al .  An  outline  of  the 
pertinent  formalism  together  with  tabulations  of  inter¬ 
mediate  steps  are  given  in  Appendix  II.  These  calculations 
require  the  radial  part  of  the  single-particle  wave 
functions.  Since  the  method  of  Glaudemans  et.  al .  did 
not  require  an  explicit  form  for  these  wave  functions 
it  was  necessary  to  assume  them.  For  simplicity 
harmonic  oscillator  wave  functions  of  the  following 
form  were  chosen. 
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The  quantity  is  a  normalization  factor,  L  is 

an  associated  Laguerre  polynomial  and  b  is  an  unknown 

parametei .  The  value  of  b  was  taken  from  the  work  of 

Arima  (Ar  58)  as  1.826  fermis .  More  specifically, 

the  normalized  U’s  for  the  2s^  and  ld^  shells  are 

1  2 

1  3 

U2s^  =  3Tr4b2(3-  2^—)  exp  (  - 1  ^-) 

•  b2  2  b2 


U 


Id 


(r) 


1  3 

4  i  2  r2 
b  —  exp 

b2 


A  summary  of  calculated  multipole  mixing  ratios, 
branching  ratios  and  lifetimes  is  presented  in  figure 
4-1  together  with  experimental  data  obtained  in  the 
present  work.  It  is  seen  that  experimental  lifetimes 
of  the  0.68  and  0.71  MeV  levels  and  branching  ratios 
of  all  states  considered  except  the  2.94  MeV  level 
agree  reasonably  well  with  the  calculated  values.  In  those 
cases  of  agreement  with  experiment,  values  calculated  using 
wave  functions  from  Glaudemans '  shell  model  are  ap¬ 
proximately  the  same  as  the  values  predicted  by  the 
Weisskopf  single-particle  estimates.  A  factor  of 
100  was  used  to  decrease  the  single-particle  estimate 
for  the  AT  =  0  transition  from  the  0.71  MeV  level  to 

\ 

the  ground  state  because  of  isospin  selection  rules. 
Branching  ratio  calculations  for  the  2.94  MeV  level 
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Figure  4-1  A  Comparison  of  Experimental  and  Calculated 

Multipole  Mixing  Ratios,  Branching  Ratios 
and  Lifetimes  for  5  Low  Lying  States  of 

30P. 
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are  in  marked  disagreement  with  experiment.  The  only 
mixing  ratio  that  agrees  reasonably  well  with  experiment 
are  those  for  the  0.71  -*  0  and  2.54  -*  0  MeV  transitions. 
For  the  2.7  2  0  MeV  transition  there  is  a  significant 

difference  between  experimental  and  calculated  mixing 
ratios.  Thus,  it  appears  that  a  shell  model  which 
assumes  an  inert  28Si  core  does  not  adequately  describe 
radiative  properties  of  the  2.72  or  2.94  MeV 

levels  of  30P. 
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CHAPTER  V  CONCLUSION 


The  following  spin  possibilities  have  been  obtained 

3  0 

for  5  previously  unassigned  excitated  states  of  P: 


Energy  (MeV) 

J 

2.839 

i — 1 

ro 

** 

(N 

3.736 

1,  (2,3) 

3.836 

1,2,3 

3.928 

1,2,3 

4.734 

1,2,3 

In  addition,  confirmation  or  consistency  with  previous 
spin  assignments  to  states  at  1.453,  1.974,  2.536,  2.722, 
3.020,  4.142,  4.501,  and  4.634  MeV  excitation  energy  has 
resulted.  However,  both  the  spin  assignment  to  and 
branching  ratios  for  the  4.921  MeV  level  are  in  serious 
disagreement  with  the  work  of  Harris  and  Hyder  (Ha  67A) . 
Moreover,  the  25%  branch  to  the  4.231  MeV  level  indicated 
by  our  results  rules  out  the  possibility  of  a  J=4  assign¬ 
ment  to  the  4.231  MeV  level  (the  assignment  made  to  the 
4.231  MeV  le'vel  by  Harris  and  Hyder)  if  our  assignment  to 
the  4.921  MeV  level  is  correct.  The  branching  ratios  from 
the  present  work  for  the  2.839,  3.736,  3.836,  3.928,  and 
4.501  MeV  levels  are  also  in  disagreement  with  those  of 
previous  investigators  (Mo  63,  Ha  67A) .  The  branching 
ratios  for  the  remaining  states  investigated  showed  general 
agreement  with  previous  results. 

A  test  of  the  shell  model  of  Gaudemans '  et  al .  (G1  64) 

in  predicting  radiative  properties  of  the  0.680,  0.709, 
2.536,  2.722,  and  2.939  MeV  levels  showed  it  to  be  in¬ 
adequate  for  describing  the  2.722,  and  2.939  Mev  levels. 

It  appears  that  a  model  which  includes  core  effects  is 
required . 


' 
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APPENDIX  I 

LEAST-SQUARES  ANALYSIS 

The  following  is  a  development  of  the  least-squares 
method  of  analysis  used  in  the  present  work. 

Equation  1-1  can  be  restated  as 


W(e)  =  l  l  U(k,n,a)  UL_  p(a)p  (cos  6)  (Al-1) 

k  a  n  =  0  1 + x 

k  even 

where 

U(k , 0 , a)  =  pk(aa) Fk(LLab)Qk 


U(k,l,a)  =  -2pk(aa)Fk(LL+lab)Qh  (Al-2) 


U(k  ,  2 , a)  =  PkCaa)Fk(L+lL+lab)Qk 


The  quantity  x  is  the  mixing  ratio ,  P(a)  is  the  population 
parameter  of  the  magnetic  substate  a  and  Pk  is  the 
Legendre  polynomial  of  order  k.  The  angular  distribution 
at  each  angle  0^  in  matrix  form  is  given  by 


,y*«| 


-15  7- 


r~ 


A 

Vv. 


W(01)\ 

W(02) 


W(0N) 


n 


n  \ 


ft-  Po(cos0n)  P ^ (cos 0 1 )  .  .  ,\  /j>U(0,n,l)-~-T.  . -.)iUC0,Ti,c)-x— 3 


1  P2(cos02)  P^  (cosOp)... 


l  P2 (coseN)  P4 (coseN) 


-I  u  2 

1+x  n 


?  ' 

1+x  : 


I 


n  n 

I U ( 4  ,n,  1)— j-.  .  .J|U(4,n,c)-^- 


1  +  X 


n 


I  ,  2 

1  +  x 


x 


X 


I  P(ax) 
P  (a2) 


P  Oc), 


(Al-3) 


or  in  more  compact  notation  by 


W  =  AX 

If  Y  is  a  column  vector  of  measured  values,  then  (Y  -  W) 
is  a  column  vector  of  residues.  In  matrix  notation, 
therefore,  the  function  x2  (equation  1-14)  can  be  written 
as 

x2  "  (Y  -  W)T  to (Y  -  W)  (Al-4) 

T 

Here  Z  denotes  the  transpose  of  Z  and  w  is  the  matrix 
of  weighting  values  which  in  tlie  present  case  is  diagonal 
with  elements  given  by 

w(I,I)  =  -9— 1 —  (Al-5) 

E-Cep 

2 

The  quantity  E  (0.)  is  the  uncertainty  in  the  number  of 
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counts  at  angle  Qj  ,  Expanding  equation  CM- 4)  gives 

X2  =  YTwY  -  yTw  -  WT<dY  +  WT0)W 
Substituting  equation  (Al-3)  in  the  above  yields 


rp  *T»  rp  rp  rp  rp 

X2  =  Y  w  Y  -  Y  uAX  -  X  A  coY  +  X  AoAX 

g  y  2 

The  set  of  equations  defined  by  — 3. —  o  can 

9  P  (  a  ) 

give 

be  shown  to 

2  [ATojAX  -  Aro)Y]  =  0 

or 

X  =  (A1ojA)  ~1A1-  03 Y 

(Al-6) 

Defining  the  normal  matrix,  N,  as 


T 

N  =  A  03A 

(Al- 7) 

gives 

- 1  T 

X  =  N  A  03Y 

(Al-6) 

It  is  a.  general  result  of  statistics  that  the 

variance- 

covariance  matrix  of  the  fitted  parameters  is 

the  inverse 

of  the  normal  matrix.  Thus  the  standard  deviation  of  the 
i^1  element  of  X  is  given  by 

0.  =  (N"bu  (Al-8) 


which  in  this  case  is  the  uncertainty  in  the  population 
parameters  P(a..  ).  After  solving  for  X  one  can  compute 
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W  for  a  given  value  of  the  mixing  ratio,  x,  from  equation 
Al-3  and,  hence,  compute  the  value  of  x2* 

In  simultaneously  fitting  two  transitions  both  of 
which  arise  from  a  given  level  decay  one  replaces  Al-3 
by  the  matrix  equation 

(Wx  +  W2)  -  (A1  +  A2)X  CA1-9) 

where  the  subscripts  1  and  2  denote  transitions  1  and 
2  respectively.  This  formulation  is  possible  since 
the  unknown  population  parameters  (contained  in  X)  are 
the  same  for  both  transitions.  Extension  to  three  or 
more  transitions  involves  replacing  both  W  and  A  of 
Al-3  by  the  sum  of  the  component  W’s  and  A's  corresponding 
to  each  transition.  The  additional  requirement (s)  of 
having  to  fit  two  or  more  distributions  will  generally 
result  in  fewer  possible  values  for  the  population 
parameters  and,  hence,  more  restricted  mixing  ratio 
solutions  together  with  better  discrimination  between 
possible  spin  values. 
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APPENDIX  II 

CALCULATION  OF  REDUCED  MATRIX  ELEMENTS  FOR 
TRANSITIONS  BETWEEN  TWO-PARTICLE  STATES 

Formulas  from  the  Rose  and  Brink  phase  consistent 
formalism  for  calculating  multipole  mixing  ratios, 
branching  ratios  and  lifetimes  for  transitions  between 
two-particle  states  are  given  in  the  following  presenta¬ 
tion.  These  formulas  are  followed  by  a  tabulation  of 
intermediate  steps  in  the  calculations  for  two -particle 
states  in  the  anc^  ^^3/2  shells.  The  symbols  are 

the  same  as  those  used  by  Rose  and  Brink  (Ro  67)  . 

The  essential  problem  in  predicting  radiative 
properties  of  transitions  between  two-particle  states 
involves  calculating  reduced  matrix  elements  of  the  form 


<  £ 


1 

1  2 


1  <  7T  > 

L 


n  1  • 

l2  2  J 


The  ket  | i  -  j>  is  the  wave  function  for  a  spin-- 


particle  with  spin- orbit  coupling  in  a  central  field 
1  -> 

(i.e.  l  +  2  =  l)  and  denotes  the  type  of  interaction 


(i.e.  electric  or  magnet 

ic.) 

The 

rp  <  IT  > 

A  LM 

are  one-particle 

operators  defined  by 

rp  0 

LM 

e 

aL 

^LM  + 

Wf-P 

(A2- 1) 

Tm 

LM 

m 

aL 

<mlm  + 

MuP 

(A2-1 ') 

—  \ 


■ 
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The  multipole  expansion  coefficients,  cxj77"  are  listed 
up  to  L  =  3  in  Table  A2-1. 


TABLE  A2-1  Multipole  Expansion  Coefficients 


L 

e 

aL 

m 

aL 

1 

ik 

k 

2 

-k2/2/3 

ik2/2/3 

3 

— 

-ik3v/2/15/3 

-  k  3  /2  / 1 5  /% 

For  convenience  Rose  and  Brink  have  replaced  the  operator 
QLm  bY  (Qpjvp  eff  (  =  eg£rLcLM)  •  Here  e  denotes  electronic 
charge,  r  is  magnitude  of  the  vector  r  which  describes 
the  particle’s  position,  g  is  zero  for  a  neutron  and  1 
for  a  proton, f  is  related  to  spherical  harmonics, 


Y 


LM 


by 


4  IT 


'LM 


c  2L  +  1) 


.  V 

h  LM 


For  a  transition  from  a  single-particle  state  | A>  to 

another  single-particle  state  | y >  by  emission  of  a  gamma 

ray  (i.e.  E  >  E  )  the  matrix  element  of  QTM  is  given  by 
A  y  i.i'i 


<xlQLMlv>  =  2S 


£ 


Sc)CEx  - 


h  k  / 


' 


■ 
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and  that  for 


is  given  by 


y>  =  eg£<r^ 


Lr 
r  (> 


LM 


y; 


The  quantity  m  denotes  the  mass  of  a  nucleon  6'  is  the  nuclear 
magneton  and  lc  is  the  photon  wave  number.  Since  the  photon 
wave  number  is  given  by 


k  = 


E  -  E 
. A  .  ]i_ 

\c 


in  this  case 


<A  IQlmI  v>  <A  I  ^Lldeff  I ll> 


However,  when  the  single- particle  energies  E,  and  E  arc 


A 


such  that  E^  <  E^  ,  one  must  change  the  sign  of  <X  |  (Q,  j  •)  c,  ,p£  |  |i> 
The  single-particle  reduced  matrix  elements  correspond¬ 
ing  to  the  operators  (QLM)ef£,  Q£M>  MLM  and  M^M  are  given 
by 

<li  2  '  ^L^eff  I  I  £2  2  ^ 2 >  =  eg£J^b12  ^A2‘2^ 


<£1  2  ^  1 1  l^iJ  I  l2  2  J  2 


t  ,  i  +  - 

<4j  \  j1l|MLIU2  l  j2>  =  2Bgp  3  2(-)'  J  2I(L-l)(i1£,0l|LX) 


l  g  gkl(L)  (a^ap-i^-  (A2 - 2  '  ) 
z  s  1  z  L+l 

1 

TfT.-lHo  , 

12 


1 

X  W(j1B]  j2i2;|  10  [  (2j  2  +  l)  (2^+n  (2il2  +  l))l2(£2H-l)L]2 


x  [  (2L+1)  (J.+  l)  ] 


1 

2 


(A2- 2” j 
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<£1  \  jlHMLH£2  1  -V  "  2  8sBI(L-l)(L-a1-a2)b12 


(A2-2,M) 


where  I (L)  is  the  radial  integral 


I(L)  =  /  U^(r)r^  ’  ^  U^,(r)dr  (A2-3) 

an  d  .. 

o-o  i--  -  _  _ 

i  .1  2  r  dl.2rr  lV12r9T  ,  n  .  1  l,Tn, 

bj2  =  1  (-)  (2j2  +  1)  (2L+1)  (j  n  22 "  2  '  ^'° 

(A2-4) 

ax  =  -  j1)(2j1  +  1)  (A2-5) 

The  quantity  gs  is  5.585  or  -3.826  for  a  proton  or  neutron, 
respectively . 

Clebs ch- Gordan  and  Racah  coupling  coefficients  are 

denoted  by  C  j  -j  j  2  £l£2^3£3^  andw(j‘1j2  £2£1;^3£3^J 
respectively.  To  conserve  parity  L+  J^+  f-,  must  be  even 
for  electric  operators  and  odd  for  magnetic  ones.  The 
U  1  s  are  radial  components  of  the  single-particle 
wave  functions.  A  tabulation  of  I(L)  values  for  harmonic 
oscillator  wave  functions  is  given  in  Table  A2 -  2 . 

Reduced  matrix  elements  for  transitions  between  two- 
particle  states  can.  be  expressed  in  terms  of  those  for 
one -particle  states  by  the  expression 


L+j  -j, 

<jajb  JT||TL(l)+TL(2)||jcjdJ'T'>  =  (-)  C  D  (2J 1  -i-l) 


1 

2 


1 


X 


1 

2 


•  ””  \ 
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TABLE  A2-2 


Summary  of  I (L)  Values  for  Harmonic  Oscillator 


Wave  Functions  as  a  Function  of  the  Parameter  b. 
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1 


1 


1 


1 


(A2-6) 


where 


(A2-7) 


=  j  r  -i-  j .  . 
'  U  D 


It  is  arbitrary  in  the  case  of  a  neutron-proton  two  body 
problem  whether  nucleon  ”1"  is  the  proton  and  ”2"  the 
neutron  or  vice  versa. 

From  two-particle  reduced  matrix  elements  calculated 
in  the  above  way  one  can  determine  reduced  matrix  elements 
for  different  two-particle  configurations.  The  predicted 
configurations  are  obtained  from  model  calculations. 

Using  the  latter  reduced  matrix  matrix  elements  one  can 
calculate  both  the  sign  and  the  magnitude  of  multipole 
mixing  ratios  from  the  definition 


(A2- 8) 


-166- 


The  symbols  L  and  tt  denote  the  lowest -order  multipolarity 
occurring  in  the  transitions  J.,  J ^  •  One  can  also 
calculate  the  speeds  for  given  transitions  and  lifetimes 
of  nuclear  states.  The  transition  speed,  t!',  for  a  given 
type  of  radiation  (i.e.  electric  or  magnetic)  with  multi¬ 
polarity  L  is  given  by 


1  4  k 


TL1'>  I  |J2>  |Z/(2L+1) 


(A2-9) 


The  total  transition  speed 
a  mixing  ratio  6  is 


t  . 
i 


for  a  given  transition  with 


—  =  (1  +  62)--t~ 
t  •  L 

1  T  • 


(A2-10J 


l 

where  L  refers  to  the  lowest  multipole,  component .  Finally, 
the  lifetime  of  a  nuclear  state  "an  is  given  by  the  sum 


(A2-11) 


l  l 

where  i  refers  to  all  transitions  that  go  from  state  "a” 
to  lower  lying  states.  Branching  ratios  are  obtained 

T 

from  the  ratios  —  . 

t  . 


l 

A  summary  of  values  for  the  various  matrix  elements 
of  equations  (A2-2)  (as  applied  to  the  2s1^2  and  ld^^ 


' 
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shells)  is  given  in  Table  A2-3.  On  this  tabulation  the 

integral  I (L)  between  states  £  and  h  is  denoted  by 

(I(L))^.  The  sign  of  the  2s^^  ld^/2  matrix  element 

of  (Q9)  r  r  has  been  cdianged  since  E9  <  E  , 

z  eu  Zsi/2  lci3/2 

Next  a  summary  of  two -particle  matrix  elements  for  both 

Ml  and  E2  multipoles  for  T  +  T5  even  is  given  in  Table 

A2-4.  A  summary  of  calculated  configurations  for  five 
30 

states  of  P  as  obtained  by  Glaudemans  e_t  a  1 .  ( G 1  6) 

is  given  in  Table  A2-5. 


f 


■I 
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TABLE  A2-3a  Proton  Single  -  Part icle  Matrix  Elements 

<£1  ^  J*  1 1  I  ^L^eff  I  I  £  2  **  ^2> 


Initial 

State 

Final 

State 

L  =  1 

L  =  2 

L  =  3 

2si 

2si 

0 

0 

0 

2 

2 

2si 

ld3 

0 

0.633  e(I(2))sd 

0 

2 

2 

ld3 

2si 

0 

0.447  e(I(2))sd 

0 

2 

2 

ld3 

ld3 

0 

-0.447  e(I(2))dd 

0 

2 

2 
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TABLE  A2-3b  Single  Particle  Matrix  Elements 

<£1  35  J  1 1  I  Ql  I  I  a2  35  V 


"  "  . — — 

Particle 

Initial 

State 

Final 

State 

L  =  1 

L  =  2 

L  =  3 

Proton 

2si 

2si 

0 

0 

0 

2 

2 

2s^ 

ld3 

0 

1.766Bk(I(2))sd 

0 

2 

2 

ld3 

2si 

0 

1.249ek(I(2))sd 

0 

2 

2 

ld3 

ld3 

0 

0 

0 

2 

7 

Neutron 

2s^ 

2si 

0 

0 

0 

2 

2 

2s-^ 

ld3 

0 

-1.210  Bk(I(2))sd 

0 

2 

2 

ld3 

2si 

0 

-0.856  Bk (I (2) ) sd 

0 

2 

2 

ld3 

ld3 

0 

0 

0 

i - 

2 

2 
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TABLE  A2-3d  Single-Particle  Matrix  Elements 


j 


2 


> 


"  - 

Particle 

Initial 

State 

F  inal 

State 

L  =  1 

L  =  2 

L  =  3 

Proton 

2si 

2si 

4.84  3 

0 

0 

2 

2 

2si 

ld3 

0 

0 

0 

2 

2 

ld3 

2s-^ 

0 

0 

0 

2 

2 

ld3 

ld3 

-2.16  6 

0 

1.416  g(I(2))dd 

2 

2 

Neutron 

2s-^ 

2si 

-3.31  6 

0 

0 

2 

2 

2si 

ld3 

0 

0 

0 

2 

2 

ld3 

2  s 

0 

0 

0 

2 

I 

ld3 

ld3 

1.48  6 

0 

-0.970  B(I(2))dd 

2 

2 

J 
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TABLE  A2-3c  Proton  Single  Particle  Matrix  Elements 


<  £  1  ^2  j  1 1  I M 


L 


£  2  35  i  2y 


Initial 

State 

Final 

•State 

L  =  1 

L  -  2 

L  =  3 

2si 

2si 

0 

0 

0 

2 

2 

2sl 

lcl3 

0 

0 

0 

2 

2 

ld„ 

2s., 

0 

0 

0 

3 

1 

2 

2 

1 a  3 

lcl3 

2 . 3246 

0 

-1.014B(I(2))dd 

2 

2 

TABLE  A2-4a  E2  Two-Particle  Matrix  Elements  for  T  +  T'  Even. 


in  the  2s j  and  ld^  Shells. 
2  2 


Final 

I  n  i  t  i  ai\ 

!  2 

Is!" 

2 

KV 

2  2 

|d2> 

2 

2 

0 

[0.0913e+0.0803ak] 

k2(I.(2))sd 

0 

1  s  d  > 

|S1  3 

2  2 

- [0.0913e+0. 08036k] 
k2(K2))sd 

0.0457ek2(I(2))dd 

[0. 0289e  +  0. 0254  3k] 
k2(l  (  2  )  )  s  d 

id?> 

1  3 

2 

0 

- [0 . 0 2 8 9 e  +  0 . 0254  3k] 
k2(I(2))sd 

-0.07301<2(I(2))d(I 

1 


-172- 


TABLE  A2-4b  Ml  Two-Particle  Matrix  Elements  for  T  +  T* 

Even  in  the  2s1  and  l.d3  Shells. 

2  2 


^N^Final 

Init  ialSvv. 

/ 

i  2 

1  Si> 

2 

| s1d3> 

2  2 

ld3>2 

2 

i  2 

|s1> 

2 

1.2443k 

0 

0 

| s1d3> 

2  2 

0 

0.4443k 

0 

'd3> 

2 

0 

- — - - 

0 

0.6003k 

i 

TABLE  A2-5  Shell  Model  Two-Particle  Configurations  of 

sn 

Predicted  P  states. 


Level  Energy  (MeV) 

# 

Amplitudes 

Expt ’ 1 

i 

Predicted 

I  si> 

KV 

14 

J 

T 

2 

2  2 

2 

1 

0 

0 

0 

+0.482 

+0.865 

+0.138 

0 

1 

0.680 

0.65 

-0.852 

+  0.0 

+0.523 

1 

0 

0.709 

0.67 

-0.876 

+0.471 

+0.100 

3 

0 

2.536 

2.71 

+0.000 

+0.000 

+1.000 

2 

0 

2.  722 

2.89 

0.000 

+1.000 

+0.000 

2 

1 

2.939 

2.98 

0.000 

-0.971 

+0.239 

■ 


. 


X 


. 
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LIST  OF  SYMBOLS 

A 

matrix  related  to  Pk's,  U(k,n,  )’s  and  x’s 

A. 

1 

corrected  photopeak  area  corresponding  to 

Ak 

coefficients  in  a  Legendre  polynomial  expan¬ 
sion  of  order  k 

e 

electronic  charge 

E 

X 

single-particle  energy  of  state  |x> 

Fk(ab) 

coefficients  in  the  expression  for  W(q) 

which  depend  only  on  the  gamma-ray  transi¬ 
tion  from  state  'a*  to  state  *bT 

Fk(LL’ab) 

coefficients  in  the  expression  for  W(6) 

related  to  the  above  F’s 

i 

imaginary  unit 

I  (L) 

integral  involving  single-particle  radial 

wave  functions 

£ 

single-particle  angular  momentum 

k 

Wave  number  (  =  2tt/x)  where  X  is  the  wave 

length  of  a  photon 

£ 

orbital  quantum  number 

L 

multipolarity  of  a  gamma-ray  transition 

Lm 

associated  Laguerre  polynomial 

M 

magnetic  quantum  number  corresponding  to  L 

|  M  |  2 

strength  of  a  transition  (gamma)  in  Weisskopf 

units 

N 

number  of  degrees  of  freedom  in  x2 

P(a) 

t  h 

population  parameter  of  a  1  magnetic  substate 

-174- 


pk(cos  0) 
Qk 

^LM’^LM 

T 

LM 

Uk(ab) 

U (k ,n , a) 

U 

n  £ 
w(e) 

X 


X 

Y(8) 


ylm 

a 


3 


Pk(a) 

Pk(aa) 


Legendre  polynomial  of  order  k 
attenuation  coefficients  for  the  gamma-ray 
detector 

electric  multipole  operators 
general  interaction  operator 
coefficients  similar  to  Fk(ab) 
coefficients  in  W(0)  related  to  F^LL'ab), 

Qk  and  Pk(ac0 

radial  part  of  harmonic  oscillator  wave  function 
angular  distribution  of  a  gamma  ray 
column  vector  whose  elements  are  the  population 
parameters 

mixing  ratio  of  a  gamma-ray  transition 

gamma-ray  yield  at  angle  0 

spherical  harmonic 

magnetic  quantum  number 

multipole  expansion  coefficient 

nuclear  magneton 

chi-square 

matrix  of  weighting  values  in  the  expression 
for  x2 

statistical  tensor  describing  the  population 
t  h. 

of  the  a  magnetic  substate  of  state  ’a’ 
coefficients  occurring  in  the  expansion  of 
Pk(a)  in  terms  of  the  population  parameters 
standard  deviation 


a 


. 
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( I  ^m) 

w(abcd ; ef ) 
<a | | H | | b> 


transition  speed  in  seconds 

the  Rose  and  Brink  notation  for  a  gamma-ray 

mixing  ratio 

Clebsch-Gordan  coefficient 
Racah  coefficient 

a  reduced  matrix  element  of  the  operator  H 
between  states  raf  and  * b ’ 
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